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ABSTRACT 


The  contribution  of  this  thesis  is  the  development  of  a  CAD  (computer  aided 
design)  tool  for  current  mode  multiple-valued  logic  (MVL)  CMOS  circuits.  It  is 
only  the  second  known  MVL  CAD  tool  and  the  first  CAD  tool  for  MVL  CMOS. 

The  tool  accepts  a  specification  of  the  function  to  be  realized  by  the  user, 
produces  a  minimal  or  near-minimal  realization  (if  such  a  realization  is  possible), 
and  produces  a  layout  of  a  programmable  logic  array  (PLA)  integrated  circuit  that 
realizes  the  given  function.  The  layout  is  i.i  MAGIC  format,  suitable  for  submission 
to  a  chip  manufacturer.  The  CAD  tool  also  allows  the  user  to  simulate  the  realized 
function  so  that  he/she  can  verify  correctness  of  design. 

The  CAD  tool  is  designed  also  to  bean  analysis  tool  for  heuristic  minimization 
algorithms.  As  part  of  this  thesis,  a  random  function  generator  and  statistics  gath¬ 
ering  package  were  developed.  In  the  present  tool,  two  heuristics  are  provided  and 
the  user  can  choose  one  or  both.  In  the  latter  case,  the  better  realization  is  output 
to  the  user.  The  CAD  tool  is  designed  to  be  flexible,  so  that  future  improvements 
can  be  made  in  the  heuristic  algorithms,  as  well  as  the  layout  generator.  T  hus, 
the  tool  can  be  used  to  accommodate  new  technologies,  for  example,  a  voltage  mode 
CMOS  PLA  rather  than  the  current  mode  CMOS  currently  implemented. 
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I.  INTRODUCTION 


A  computer-aided  design  (CAD)  tool  has  been  created  at  the  Naval  Postgrad¬ 
uate  School  to  design  multiple- valued  logic  programmable  logic  arrays  (MVL-PLA). 
The  CAD  tool  accepts  a  function  specification  and  produces  a  circuit  layout.  Small 
MVL-PLA's  can  be  designed  by  hand.  However,  for  most  practical  applications, 
hand  design  is  too  time  consuming  and  prone  to  errors. 

The  CAD  tool  performs  the  following: 

•  Accepts  an  input  function  specification,  specifically  an  algebraic  expression. 

•  Minimizes  the  expression  producing  a  realization  that  occupies  the  least 
amount  of  space. 

•  Produces  a  layout  which  can  be  used  to  generate  the  MVL-PLA  chip. 

The  CAD  tool  is  designed  toaccommodate future  improvements  in: 

•  The  minimization  program  and 

•  MVL-PLA  layout. 

It  is  expected  that  as  knowledge  is  gained  about  these  aspects  of  MVL-PLA 
design,  there  will  be  new  subroutines  that  will  substitute  for  these  parts  of  tlm  tool. 
The  whole  program  has  been  written  so  that  such  changes  can  be  easily  made. 

This  CAD  tool,  which  applies  to  current-mode  CMOS  technology,  is  only  the 
second  known  tool  for  multiple- valued  logic  circuit  design.  A  CAD  tool  exists  for 
the  design  of  multiple-valued  CCD  PLA  design  [Ref.  1). 
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A  BACKGROUND 

The  promise  of  current-mode  CMOS  circuits  for  compact  realizations  has  gen¬ 
erated  much  interest.  Japanese  researchers  have  developed  a  binary  multiplier  in¬ 
tegrated  circuit  using  bi-directional  current-mode  multiple-valued  CMOS  [Ref.  2]. 
The  circuit  has  speed  almost  equal  to  the  fastest  binary  multiplier  [Ref.  3]  but 
requires  half  the  chip  area  and  dissipates  half  the  power.  Researchers  in  Holland, 
France  and  the  United  States  have  also  fabricated  current  mode  multiple- valued 
CMOS  circuits. 

Since  this  technology  is  new,  no  design  tool  has  been  developed  for  it.  Design 
tools  are  absolutely  needed  to  produce  practical  circuits  since  modern  VLSI  circuits 
are  so  complicated.  It  is  the  lack  of  such  tools  that  has  motivated  the  research 


described  in  this  thesis. 


II.  PROGRAM  DEVELOPMENT 

The  package  developed  in  this  thesis  is  a  comprehensive  program  written  in 
the  “C"  language  to  run  on  the  VAX-1 1/785  computer  at  the  Naval  Postgraduate 
School.  As  “C”  is  highly  structured,  the  program  is  easy  to  maintain  and  easy  to 
develop  into  a  large  program  package.  Also,  use  of  “C”  assures  compatibility  with 
presently  available  binary  VLSI  design  tools  such  as  the  Berkeley  CAD  tools  [Ref. 

4]- 

A.  APPROACH  TO  THE  PROBLEM 

A  top-down  programming  technique  was  used.  First,  the  tasks  were  defined 
which  the  program  needed  to  perform.  Then,  each  task  was  developed  through  a 
series  of  steps  to  identify  very  specific  routines  required  to  perform  each  task. 

The  major  tasks  are  as  follows: 

•  Data  Entry 

•  Logic  Minimization 

•  Logic  Simulation 

•  Layout  Generation 

1.  Data  Entry 

—  Allows  the  user  to  enter  input  expression  from  an  input  data  file. 

-  Partitions  input  expression  into  identifiable  parts,  called  tokens. 

-  Examines  each  token  foi  errors. 

-  Flags  errors  when  they  occur. 
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2.  Logic  Minimization 

-  Three  minimization  algorithms  are  used. 

*  Pomper  and  Armstrong  [Ref.  5]. 

*  Dueck  and  Miller  [Ref.  6j. 

*  Gold  [Ref.  7], 

—  All  of  them  are  heuristic  and  based  on  the  direct  cover  approach. 

—  Minimization  produces  a  minimal  or  near  minimal  solution,  thus  reducing 
the  chip  area  needed. 

This  routine  will  be  discussed  in  more  detail  in  the  following  section. 

3.  Logic  Simulation 

-  Verifies  the  minimal  solution  from  the  Logic  Minimization  routine  by 
applying  test  vectors  given  by  the  user  to  it. 

4.  Layout  Generation 

-  Uses  a  minimal  solution  from  the  Logic  Minimization  routine  to  produce 
a  layout  which  will  be  used  for  MVL-PLA  chips. 

B.  ORGANIZATION  OF  THE  PROGRAM 

The  basic  organization  of  the  program  is  shown  in  Figure  2.1.  The  solid  lines 
show  flow  of  control.  The  user  supplies  input  expressions  via  the  data  entry  routine. 
These  are  applied  to  the  logic  minimization  routine  for  obtaining  the  minimized 
expression.  The  minimized  expression  is  then  passed  to  the  logic  simulation  routine 
for  verification.  Then  it  is  used  for  layout  generation  for  MVL-PLA  chips. 
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A  CAD  TOOL 


Figure  2.1:  Program  Organization 

PLA  design  using  this  program  resembles  program  development  under  an  op¬ 
erating  system.  An  editor  is  used  to  develop  the  source  file  (input  expression).  This 
is  compiled  by  the  data  entry  routine.  If  it  is  free  of  syntax  errors,  it  is  assembled 
(minimized)  and  run  (a  layout  is  produced). 


III.  PROGRAM  DESCRIPTION 


The  CAD  tool  is  best  described  by  considering  a  specific  example.  The  user 
starts  with  a  logic  diagram  (or  Karnaugh  map)  of  the  input  expression  to  be  realized. 
An  example  of  the  logic  diagram  is  shown  in  Figure  3.1.  It  shows  a  4- valued  2 
variable  expression.  Without  a  CAD  program,  an  inexperienced  user  might  spend 
half  an  hour  to  find  a  minimal  solution  of  the  expression.  What  is  worse,  when 
the  expression  is  complicated  because  of  a  large  radix  and/or  a  large  number  of 
variables,  it  may  be  impossible  not  only  to  find  the  minimal  solution,  but  also  to 
visualize  it. 

Since  the  CAD  program  is  designed  to  find  a  minimal  solution  of  the  expres¬ 
sion,  the  user  is  not  required  to  minimize  the  expression.  All  the  user  has  to  do  is 
specify  the  expression  correctly.  Figure  3.2  shows  one  way  the  function  of  Figure’ 
3.1  can  be  expressed  for  the  CAD  tool.  Here  each  circle  represents  a  product  term. 
Equation  3.1  shows  an  algebraic  equivalent  of  the  function  of  Figure  3.2. 


4:2:  (3.1) 

+2  *  AT (1,2)  *  A’2(0,0) 

+3  *  AT(3,3)  *  A’2(l,  1) 

+2  *  AT  (0,0)  *  ,Y2(1,2) 

+  1  *  AT (0,0)*  *2(2,3) 

+2  *  AT(1, 1)  *  A’2(l,  1 ) 

+2*  AT ( 1 , 1 )  +  A'2(2,3) 
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+  1  *  Yl(l,  1)  *  A'2(2,2) 

+2*  A'1(2,3)*.Y2(1,1) 

+  1  *  A' 1(2, 2)  *  A'2(l,2) 

+  1  *  A'1(3,3)*.Y2(2,2) 

+3  *  A' 1(2, 2)  *  A’2(3, 3) 

Equation  3.1  is  represented  according  to  the  standard  format  of  the  expression 
shown  in  Table  3.1.  Any  multiple- valued  logic  function  f(X  1,  A’2,  A’3,  •  •  • ,  A  n)  can 
be  represented  by  the  following  format: 

/  =  radix  :  number  of  variables  :  (3.2) 

+Ci  *  X 1  (a  1, 61)  *  A’2(«2, 62)  *  •  •  •  *  Xn(an,  bn) 

+ e2  *  XT(cl,rfl)  *  A"2(c2, d‘2)  *  ■  ■  ■  *  Xn(cn,dn ) 

+c„  *  Xl(fl,gl)  *  X2(f2,g2)  *■■■*  Xn{fn,gn) 

where  Cj,  c2,  •  •  •  and  cn  denote  a  logic  value  between  ‘I’  and  ‘(radix  -1 )',  the  symbol 
‘*’  denotes  the  minimum  operator,  Xn(an,  bn)  denotes  a  literal  function  of  the  input 
variable  ‘Xu’  with  lower  bound  ‘an’  and  upper  bound  ‘bn'.  Specifically,  Xn(an,bn) 
is  (radix  -1)  if  an  <  An  <  bn,  and  0  otherwise.  The  symbol  ‘-f’  represents  the 
truncated  arithmatic  sum  (truncated  to  (radix  -1)  if  the  actual  sum  exceeds  (radix 

-  D). 

According  to  Equation  3.1,  the  expression  has  radix  4  and  2  variables  and  a 
li  *  of  terms  (11  terms  in  this  expression).  A  term  is  composed  of  a  coefficient  and 
a  list  of  variables  (2  variables  for  each  term  in  this  case).  The  expression  shown  in 
Equation  3.1  is  entered  into  the  CAD  system  via  the  data  file. 
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ified  by  the  User 


Expression 

fladlx  :  Number  of  Variables  :  A  Usf  of  lerms 

Term 

Coefficient  *  A  List  of  Variables 

Variable 

X  (  Variable  ID  )  (  Lower  Bound.  Upper  Bound  ) 

TABLE  3.1:  Standard  Format  of  an  Input  Expression  for  a  CAD  Tool 
A.  DATA  ENTRY 

Figure  3.3  shows  the  flowchart  of  the  Data  Entry  routine.  The  input  expression 
from  the  input  data  file  is  applied  to  this  routine  first.  It  examines  all  the  expressions 
in  the  data  file.  There  is  a  special  function,  called  the  parser/scanner  in  this  routine. 
It  partitions  the  expression  into  identifiable  parts,  called  tokens.  Then,  it  passes  each 
token  to  the  next  three  consecutive  error  checking  subroutines  -  the  syntax  error 
checker,  the  grammar  error  checker,  and  the  range  error  checker. 

If  an  error  is  found,  the  corresponding  error  message  is  presented  to  the  user 
and  the  program  execution  stops.  The  user  needs  to  go  back  to  the  input  data 
file  and  correct  it.  If  all  the  expressions  are  correct,  all  are  stored  in  the  allocated 
memory  space.  Figure  3.4  shows  the  data  structure  of  the  input  expression.  It.  is 
developed  to  accommodate  the  input  expression. 


10 


Input  Data  File 


INPUT 

EXPRESSION 


Figure  3.3:  Flow  Chart  of  the  Data  Entry  Routine 


2 


According  to  the  Figure  3.4,  three  data  structures  form  a  complete  input 
expression.  They  are  the  Expression,  Term  and  Bounds  data  structures.  The  radix, 
number  of  variables  and  number  of  terms  in  the  input  expression  are  saved  in  the 
Expression  structure. 

The  size  of  the  array  for  the  term  lists  is  determined  by  the  number  of  terms  of 
the  expression.  Also,  the  size  of  the  bounds  array  is  equal  to  the  number  of  variables 
used  in  the  expression.  Therefore,  all  the  components  of  the  input  expression  are 
saved  into  this  data  structure. 

Once  this  data  structure  for  the  input  expression  is  formed,  it  will  be  used 
throughout  the  CAD  program. 

B.  LOGIC  MINIMIZATION 

The  data  structure  of  the  input  expression  is  formed  in  the  Data  Entry  routine. 
This  is  then  used  by  Logic  Minimization  routine  to  minimize  the  expression.  Figure 
3.5  shows  the  flowchart  of  the  Logic  Minimization  routine. 

First,  a  working  copy  (f)  of  the  original  input  expression  is  made.  If  f  is  fully 
covered  (made  up  of  only  zeroes  and  don’t  care  terms  (radix)),  the  program  stops. 
The  minimization  process  is  accomplished.  Otherwise,  the  program  finds  the  most 
isolated  minterm  first,  then  picks  the  best  implicant  (or  product  term)  from  among 
the  possible  implicants  which  cover  the  chosen  most  isolated  minterm. 

The  best  implicant  chosen  is  saved  as  a  part  of  the  minimal  solution  of  the 
input  expression.  Then  it  is  taken  out  from  the  working  copy  (f).  This  process  is 
repeated  until  the  resulting  expression  has  no  more  minterms. 

As  mentioned  earlier,  three  heuristic  algorithms  are  used  in  this  thesis.  As 
shown  in  Reference  7,  neither  that  of  the  Dueck  and  Miller  [Ref.  6]  nor  that  of  and 
Pomper  and  Armstrong  [Ref.  5]  is  consistently  better  than  the  other.  That  is,  there 
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Figure  3.5:  Flowchart  of  Logic  Minimization  Routine 


are  functions  where  the  first  does  better  than  the  second  and  vice  versa.  Thus,  the 
technique  of  Gold  is  introduced  in  Reference  7.  This  is  simply  the  application  of 
both  algorithms  followed  by  a  choice  of  the  best  realization.  Since  the  CAD  tool  is 
also  used  to  analyze  minimization  algorithms,  implementation  of  all  heuristics  is  an 
integral  part  01  the  program.  All  of  them  work  in  an  identical  manner  except  in  two 
ways.  One  is  the  way  a  minterm  is  chosen,  and  the  other  is  the  way  an  implicant  is 
chosen. 

•  Pomper  and  Armstrong  [Ref.  5] 

The  minterm  is  chosen  randomly.  However,  the  implicant  is  chosen  as  the 
one  which  drives  the  most  minterins  to  zero  when  subtracted.  If  more  than  one 
implicant  is  available,  the  largest  is  chosen.  If  there  are  more  than  one  largest,  the 
first  generated  is  chosen. 

•  Dueck  and  Miller  [Ref.  6] 

In  this  algorithm,  an  isolation  factor  (IF)  is  calculated  for  each  minterm  with 
the  smallest  value  beginning  with  all  1  minterms,  2  minterins,  etc.  The  minterm 
with  the  largest  IF  is  chosen  and  if  there  are  more  than  one,  the  first  generated  is 
chosen.  IF  is  calculated  as  the  inverse  of  the  clustering  factor  where  the  clustering 
factor  is  as  follows: 

Cl'a  =  DEAa  x  ( radix  —  1)  +  EAa  (3.3) 

where  DEA  is  the  number  of  variables  (directions)  in  which  a  minterm  (a)  can 
be  combined  with  a  non-zero  number  of  minterms  (called  direction  of  expandable 
adjacency).  EA  is  the  number  of  adjacent  minterms  with  which  a  minterm  («)  can 
be  combined  in  an  interval  literal  (called  expandable  adjacency).  CF  provides  a 
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measure  of  the  degree  of  which  a  specific  minterm  can  combine  with  other  minterms 
in  the  expression. 

The  minterm  with  the  smallest  CF  is  chosen  as  the  most  isolated  minterm. 
All  implicants  that  cover  this  minterm  are  then  generated.  A  measure  called  the 
relative  break  count  (RBC)  is  calculated  for  each.  This  provides  a  measure  of  the 
degree  to  which  the  expression  is  simplified  if  the  implicant  under  consideration  is 
chosen.  The  implicant  with  the  smallest  RBC  is  chosen  as  the  best  implicant.  If 
there  are  more  than  one  best  implicant,  the  first  generated  is  chosen. 

•  Gold  [Ref.  7] 

Gold  is  a  heuristic  in  which  the  algorithms  of  Pornper  and  Armstrong  and 
Dueck  and  Miller  are  applied,  and  the  best  realization  (the  one  with  the  smallest 
number  of  product  terms)  is  chosen.  It  was  inspired  by  the  observation  that  these 
algorithms  display  a  diversity  in  realizations  and  one  algorithm  is  not  consistently 
better  than  the  others  over  all  expressions.  However,  certain  classes  of  expressions 
do  show  that  one  algorithm  does  perform  consistently  better.  With  Gold,  the  com¬ 
bination  of  algorithms  takes  advantage  of  the  best  features  of  each. 

Table  3.2  shows  the  comparison  of  three  heuristic  algorithms.  When  logic 
minimization  is  accomplished,  the  minimal  solution  of  the  original  input  expression 
is  as  shown  in  Figure  3.6  and  Equation  3.4.  Note  that  the  number  of  product  terms 
in  the  input  expression  in  Equation  3.1  was  11  (eleven).  Now  the  minimal  solution 
(Equation  3.4)  has  8  (eight)  product  terms  in  it.  There  is  about  25%  reduction  in 
the  number  of  product  terms. 
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A1  gor  i  thin 

Choice  of 
Mi nter m 

Choice  of 
lmpl i caul 

i.  Pomper  and 
Armstrong 

Random 

Drives  Most 

Min  tot  ms  (o  0 

2.  Dueck  and 

Miller 

Largest  IT 
( Isolation 
f  actor ) 

9mal lest  RBC 

(relative  break 
count  1 

3.  Gold 

Best  of  i  and  2 

TABLE  3.2:  Comparison  of  Minimization  Algorithms 


+  1  *  AT (0,2)  *  *2(2,3) 
+  1  **1(3,3)  **2(0,2) 
+  1  *  *1(1,1)  *  *2(3,3) 
+  1  *  *1(2,3)  *  *2(1,1) 
4-2  *  *1(2,  2)  *  *2(3, 3) 
+2  *  *1(1,3)  *  *2(0,0) 
+2  *  *1(0  2)  *  *2(1, 1) 
+  2  *  *1(0, 1)  *  *2(2, 2) 
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C.  LOGIC  SIMULATION 


This  routine  is  used  to  verify  the  minimal  solution  from  the  logic  minimization 
routine.  The  user  can  specify  certain  specific  values  or  range  of  values  for  each 
variable  in  the  input  expression  (called  test  vector).  The  logic  simulation  routine 
evaluates  the  minimized  function  value  for  each  coordinate  from  the  test  vector. 
The  user  can  do  the  verification  test  by  comparing  the  original  expression  with  the 
output  of  the  Logic  Simulation  routine. 

D.  LAYOUT  GENERATION 

The  output  of  the  Logic  Minimization  routine  is  used  to  produce  a  layout  file 
(called  a  magic  file)  for  MVL-PLA  chips.  This  routine  is  discussed  in  more  detail 
in  Reference  8. 
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IV.  CONCLUSION 


This  thesis  presents  a  CAD  tool  for  design  of  multiple- valued  current -mode 
programmable  logic  arrays.  It  was  designed  to  be  user  friendly,  requiring  little 
knowledge  of  programming  techniques.  Appendix  A  contains  a  user  manual  written 
for  the  first  time  user.  In  the  interest  of  clarity,  details  of  the  program  have  been 
omitted.  However,  for  the  interested  user,  the  program  has  been  extensively  docu¬ 
mented.  This  documentation  and  the  program  are  contained  in  Appendix  B.  Three 
main  conclusions  have  resulted  from  this  work. 

•  Conclusion  1  -  Design/ Analysis  Tool 

The  CAD  tool  can  be  used  for  two  purposes.  One  is  as  a  design  tool  and 
the  other  as  an  analysis  tool.  It  is  used  for  designing  practical  current-mode  MVL 
CMOS  circuits  and  also  used  for  analyzing  and  comparing  the  different  minimization 
algorithms.  In  fact,  it  is  designed  to  get  the  statistics  for  each  minimization  program, 
which  is  very  useful  to  obtain  insight  into  logic  minimization  algorithms. 

•  Conclusion  2  -  Further  Improvement  in  CAD 

Further  research  and  improvements  are  needed  in  two  specific  areas.  These 
are  1)  Logic  Minimization  and  2)  Layout  Generation.  Still  there  are  cases  which 
the  heuristics  give  poorer  realizations  than  supplied  the  user.  Programs  are  desired 
which  can  find  the  minimal  solution  closest  to  the  optimal  solution.  In  MVL-PLA 
layout,  there  is  some  waste  of  the  chip  area.  The  further  optimized  layout  configu¬ 
ration  should  be  considered. 
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•  Conclusion  3  -  Execution  Speed  Improvement 

The  Profiler  Utility  shows  how  many  times  each  routine  is  called  and  the 
percentage  of  time  spent  in  executing  that  routine.  Below  is  an  example  output  from 
this  routine.  According  to  the  profiler  output  shown  in  Table  4.1,  there  are  several 
subroutines  which  are  called  most  frequently  throughout  the  program  execution  such 
as  _eval,  .vcopy  ,  and  _next_coord.  Typical  runs  show  these  three  routines  as  the 
most  time  consuming.  If  those  subroutines  are  carefully  examined  and  improved, 
the  total  execution  time  will  be  reduced  significantly. 
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%time 

cumsecs 

I  call 

ms/ call 

name 

92.0 

141.83 

863184 

0.16 

eval 

1.3 

143.79 

063184 

0.00 

_vcopy 

2.1 

147.00 

443878 

0.01 

next  coord 

0.0 

147.13 

6133 

0.00 

next  impl leant 

0.0 

147.16 

1789 

0.01 

copy  implicant 

0.0 

147.17 

391 

0.03 

yylook 

0.0 

147.18 

358 

0.03 

_yylex 

0.0 

147.10 

358 

0.00 

next  token 

0.0 

147.19 

356 

0.03 

streat 

0.0 

147.20 

244 

0.04 

stremp 

0.0 

147.20 

244 

0.00 

match 

0.0 

147.20 

170 

0.00 

strlen 

0.0 

147.20 

164 

0.00 

_ filbuf 

0.0 

147.20 

162 

0.00 

doscan 

0.0 

147.20 

162 

0.00 

sscanf 

0.0 

147.20 

158 

0.00 

_malloc 

0.0 

147.20 

140 

0.00 

free 

0.0 

147.20 

112 

0.00 

alloc  bound 

0.0 

147.20 

96 

0.00 

limit 

0.0 

147.20 

80 

0.00 

alloc  term 

0.0 

147.20 

80 

0.00 

_var_id 

0.0 

147.20 

80 

0.00 

var  list 

0.0 

147.20 

80 

0.00 

variable 

0.0 

147.20 

77 

0.00 

realloc 

0.0 

147.20 

38 

0.00 

bcopy 

0.0 

147.20 

32 

0.00 

gen  bounds 

0.0 

147.20 

32 

0.00 

init  implicant 

0.0 

147.20 

32 

0.00 

subtract_lmpl leant 

0.0 

147.20 

32 

0.00 

term 

0.0 

147.20 

32 

0.00 

term  list 

TABLE  4.1:  The  Output  of  the  Profiler 


APPENDIX  A 

Getting  Started 

A  few  examples  are  included  In  this  appendix  A  in  order  to  demonstrate 
the  use  of  this  CAD  tool, the  appearance  of  a  terminal  session,  and  the 
capability  of  the  utility  functions. 


1.  DESIGN  PROCEDURES 

(a)  LOGGING  IN 

-  LOG  IN  ON  A  VAX  TERMINAL. 

-  CHANGE  TO  SUBDIRECTORY  ’mag’  IN  YOUR  WORKING  DIRECTORY. 

*  cd  mag 

(b)  CREATING  THE  INPUT  DATA  FILE 

-  BY  USING  ' vi ’ EDITOR,  MAKE  AN  INPUT  DATA  FILE. 

*  vi  tt2.invl 


(1)  WHEN  THE  USER  HAS  A  SINGLE  INPUT  EXPRSSION 


4  j  2r 

+2*X1(1,3)*X2(1,1) 

+1*X1(1,3)*X2(1,3) 

+2‘X1(1,2)‘X2(0,2) 

+1‘X1(1,2)*X2(3,3) 

+2*X1(0,2)‘X2(0,2) 

+2*X1(3,3)*X2(2,3) 

+3*X1(0,2)*X2(0,0) 

+1*X1(1,2)*X2(0,2) 

+1*X1(0,3)*X2(1,2) 

+1*X1(1,2)*X2(1,2) 

+1*X1(1,2)*X2(0,1); 

(  tt2.mvl  is  the  name  of  input  data  file.  In  this  example  there  is 
only  one  expression.  It  is  a  4-valued  2-variable  expression  with 
11  terms  in  it.  A  /(semicolon)  terminates  the  expression.) 

(2)  WHEN  THE  USER  HAS  MULTIPLE  EXPRESSIONS, THE  INPUT  FILE  MAY  LOOK  AS 

FOLLOWS.  THERE  ARE  THREE  EXPRESSIONS  AND  ALL  OF  THEM  ARE  4-VALUED 

2-VARIABLE  LOGIC  EXPRESSIONS. 

4:2: 

+2*X1 (3,3)*X2(1,3) 

+3‘X1(1,3)*X2(1,2) 

+1*X1(1,3)*X2(1,1) 

+1*X1(1,3)*X2(1,3) 

+2*X1(1,2)*X2(0, 2) / 

4:2: 

+1*X1(1,2)*X2(3,3) 

+2*X1(0,2)*X2(0,2) 

+2*X1(3,3)*X2(2,3);  , 

4:2: 

+1*X1(1,2)*X2(1,2) 
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+1‘X1(1,2)*X2(0,1) 
+2*X1(0,0)*X2(1, 1) 
+3*X1(0,1)*X2(1,3); 

(input  data  file;  tt3.mvl  ) 


(3)  WHEN  THE  USER  APPLIES  TEST  VECTORS  TO  EACH  VARIABLE  FOR  VERIFICATION, 

THE  INPUT  DATA  FILE  LOOKS  AS  FOLLOWS. 

4:2: 

+2*X1(1,3)*X2(1,1) 

+1*X1(1 ,3)*X2(1,3) 

+2*X1(1,2)*X2(0,2) 

+  1»X1 (1, 2)*X2(3, 3  ) 

+2*X1(0, 2)*X2(0,2) 

+2*X1(3,3)*X2(2,3) 

+3*X1(0,2)*X2(0,0) 

+1*X1(1,2)*X2(0,2) 

+1*X1(0,3)*X2(1,2) 

+1*X1(1,2)*X2(1,2) 

+  1*XI( 1, 2)*X2(0, 1)  ; 

XI (0,3) X2 (0,3); 

(test  vector  for  the  it  out  expression  above.  Be  sure  to  put  at  the 

end.  No  need  to  put  between  variables.  In  this  case  test  vector 
has  the  range  0  through  3  for  both  variables  XI  and  X2.  The  user  can 
specify  any  arbitrary  test  vector  to  each  variable,  for  instance, 
X1(0,1)X2(2,2). ) 

(  input  data  file;  tt4.mvl  ) 

(4)  WHEN  THERE  ARE  DON'T  £ARE  TERMS  IN  THE  EXPRESS ION, USE  ' ? ' (QUESTION  MARK) 
FOR  THE  COEFFICIENT, AS  THE  FOLLOWING  EXAMPLE. 

4:2: 

+2*X1 ( 1, 3 ) *X2(1, 1) 

+1*X1(1,3)*X2(1,3) 

+2‘X1(1,2)*X2(0,2) 

+?*X1(1,2)*X2(3,3) 

+2*X1(0,2)*X2(0,2) 

+2*X1(3,3)*X2(2,3) 

+3*X1 (0,2)*X2(0,0) 

+1*X1(1,2)*X2(0,2) 

+1*X1 (0,3)*X2(1,2) 

+1 *X1 ( 1 , 2) *X2(1 , 2) 

+1*X1(1,2)*X2(0,1); 

XI ( 0, 3 ) X2( 0 , 3 ) ; 

(  input  data  file;  tt5.mvl  ) 

(  Don’t  care  term  is  4th  term  from  above  with  its  coefficient  '?') 


(C)  HOW  TO  PERFORM  LOGIC  MINIMIZATION  AND  SELECT  OUTPUT  OPTIONS  TO  GET  THE 
DESIRED  RESULTS. 

GETTING  HELP: THERE  IS  A  HELP  ROUTINE  FOR  THE  USERS  WHO  ARE  NOT  FAMILIAR 
WITH  ALL  OPTIONS  AVAILABLE  IN  THIS  CAD. 

TYPE  THE  FOLLOWING  COMMAND  FOR  HELP 
%  mvl  —  » 
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(The  following  HELP  messages  appears  on  the  screen . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 

usage:  mvl  (-eimrsvqDPGJ  (-0.0]  (-oFileJ  [-xF'ile]  (file  ...] 

-e  -  Print  the  original  expression,  as  parsed 

-i  -  Print  each  implicant  found 

-m  -  Print  the  Karnaugh  map  of  the  original  and  each 

successive  expression 
WARNING:  impractical  for  nvar  >  3 
-v  -  Verify  the  minimal  solution 

-s  -  Print  statistics  on  selected  heuristics 

-q  -  Quiet,  don't  print  final  results 

-0.0  -  Build  source  expressions  for  all  solutions 

whose  ratio  to  the  original  number  of  terms 
exceeds  this  number  and  output  them  to  a  file 
(default  is  "x.mvl") 

-D  -  Execute  the  Dueck  &  Miller  heuristic  (default) 

-P  -  Execute  the  Pomper  s  Armstrong  heuristic 

-G  -  Execute  the  Gold  heuristic 

-xFile  -  Output  source  expressions  built  using  -0.0  to  "File” 
instead  of  "x.mvl" 

-oFile  -  Ouput  formatted  solutions  to  "File” 

(Note;  In  the  above  HELP  messages  "print"  means  "output  to  CRT".) 


(EXAMPLE  1)  USING  POMPER  AND  ARMSTRONG  ALGORITHM  ON  tt2.mvl,  OUTPUT  ORIGINAL 

EXPRESSION  AS  PARKED  AND  FINAL  RESULT. 

(command  line  Is  as  follows . ) 

ece: /work/mag 
%  mvl  -Pe  tt2.mvl 

(  the  following  appears  on  the  screen . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 


Expression: 

radix:  4 

nvars:  2 

nterms:  11 
coe  f  f :  2 
XI  ( 1 , 3 ) 
X2( 1,1) 
coe  f  f :  1 
Xl( 1 , 3  ) 
X2( 1,3) 
coe  f  f :  2 
XI ( 1 , 2 ) 
X2(0, 2) 
coeff:  1 
XI ( 1 , 2 ) 
X2 ( 3, 3 ) 
coeff:  2 
X1(0, 2) 
X2(0, 2) 
coe  f  f :  2 
XI  ( 3 , 3 ) 
X2(2,3) 
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coe  f  f :  J 
Xl( 0 , 2 ) 

X2 ( 0 , 0 ) 
coeff:  1 
XI  ( 1 , 2 ) 

X2( 0, 2 ) 
coe  f  f :  1 
Xl( 0 , 3 ) 

X2( 1,2) 
coe  f  f :  1 
XI  ( 1 , 2 ) 

X2  ( 1 , 2 ) 
coe  f  f :  1 
XI  ( 1 , 2 ) 

X2(0, 1 ) 

PRA;  3/11  linplicants  -  0.27 

(  the  last  line  shows  that  the  POMTER  AND  ARMSTRONG  heuristic  produced  a 
realization  of  3  product  terms  starting  with  11  product  terms,  for  a 
reduction  of  27%) 


(EXAMPLE  2)  USING  DUECK  AND  MILLER  ALGORITHM  ON  tt2.mvl,  OUTPUT 

EACH  IMPLICANT  FOUND  AND  FINAL  RESULT. 

(command  line  is  as  follows  ...  ) 
ece:/wcrk/mag 
%  mvl  -Di  tt2.mvl 

#  . 

(  the  following  appears  on  the  screen  . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 


DRM 

M1M: 

2»X1( 

1,  1 ) *X2 ( 

3,  3) 

Imp: 

2*X1( 

1,  2 ) *X2( 

3,  3) 

DRM 

MIM: 

3  *X1  ( 

3,  3 ) *  X2( 

3,  3) 

Imp: 

3*X1( 

3,  3 ) *X2( 

1,  3) 

DRM 

MIM: 

3*X1  ( 

2,  2 ) *X2( 

0,  0) 

Imp: 

3*X1( 

0,  2)*X2( 

0,  2) 

DSM: 

3/11 

impllcants 

-  0.27 

(EXAMPLE  3)  USING  GOLD  ALGORITHM,  OUTPUT  THE  KARNAUGH  MAP  OF  THE  ORIGINAL  AND 

EACH  SUCCESSIVE  EXPRESSION  AND  FINAL  RESULT. 

(command  line  is  as  follows  . . . ) 
ece : /work/mag 
%  mvl  -Gm  tt2.mvl 

(  the  following  appears  on  the  screen . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 

Orlg  map  (DsM)i 
3.  3.  3.  0 
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3. 

3. 

3. 

3. 

3. 

3. 

3. 

3. 

0 

2 

2 

3. 

3. 

3. 

3. 

0 

3. 

3. 

3. 

3. 

3. 

3. 

3. 

3. 

0 

0 

0 

3. 

3. 

3. 

3. 

0 

3. 

3. 

3. 

4. 

3. 

3. 

3. 

4  . 

0 

0 

0 

4  . 

4. 

4. 

4. 

0 

4  . 

4  . 

4  . 

4  . 

4. 

4  . 

4. 

4  . 

0 

0 

0 

4. 

Orlg 

map  ( PSA ) t 

3. 

3. 

3. 

0 

3. 

3. 

3. 

3. 

3. 

3. 

‘  3. 

3. 

0 

2 

2 

3. 

4  . 

4. 

4. 

0 

4. 

4  . 

4. 

3. 

4. 

4  . 

4. 

3. 

0 

2 

2 

3. 

4. 

4. 

4  . 

0 

4. 

4. 

4. 

4. 

4  % 

4. 

4. 

4  . 

4  . 

0 

2 

2 

4. 

4. 

4. 

4  . 

0 

4. 

4. 

4. 

4. 

4. 

4. 

4. 

4. 

0 

0 

0 

4. 

Gold(- 

3/11 

implicants  -  0.27 

(Note;  . 

(dot) 

after  a 

number  denotes  a  highest  logic  value  or  don't 

care  generated  In  the  course  of  a  minimization.  Don't  cares  in  the 
original  function  have  no  dot. ) 


(EXAMPLE  4)  USING  POMPER  AND  ARMSTRONG  ON  tt2.mvl,  OUTPUT  THE  ORIGINAL 
EXPRESSION  AND  K-MAP  OF  THE  ORIGINAL  EXPERSSION  AND  EACH 
SUCCESSIVE  EXPRESSION  WITH  FINAL  RESULT. 

(  commnad  line  is  as  follows  . ) 

ece; /work /mag 
%  mvl  -Pei  tt2.mvl 

(  the  following  appears  on  the  screen . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 
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Expression : 


radix:  4 

nvars:  2 

n  ter  ins:  11 
coe  f  f :  2 
XI ( 1 , 3 ) 
X2( 1,1) 
coeff:  1 
XI ( 1 , 3 ) 

X  2  ( 1 , 3  ) 
coe  f  f :  2 

XI ( 1 , 2 ) 
X2 ( 0 , 2 ) 
coe  f  f :  1 
XI ( 1 , 2 ) 
X2 ( 3 , 3 ) 
coe  f  f :  2 
X1<0, 2) 
X2( 0 , 2 ) 
coe  f  f :  2 

XI  (  3 , 3  ) 
X2  (  2 , 3  ) 

coe  f  f :  3 
XI ( 0 , 2 ) 
X2 ( 0, 0 ) 
coe  f  f :  1 
XI ( 1 , 2 ) 
X2 ( 0 , 2 ) 
coe  f  f :  1 
XI ( 0, 3} 
X2 ( 1 , 2 ) 
coeff:  1 
X 1  ( 1 , 2  ) 
X2( 1 , 2 ) 
coeff:  1 
XI ( 1 , 2) 
X2 ( 0, 1 ) 


PSA 

MIM: 

3  *X1  ( 

1, 

1 )  *X2 ( 

0,  0) 

Imp: 

3*X1( 

0, 

2 )  *X2  ( 

0,  2) 

PSA 

MIM: 

3*X1  ( 

3, 

3)‘X2( 

3,  3) 

Imp: 

3  *X1  ( 

3, 

3 )  *X2( 

1/  3) 

PSA 

MIM: 

2‘X1  ( 

1, 

1 )  *X2 ( 

3,  3) 

Imp: 

2»X1( 

1, 

2)*X2( 

0,  3) 

PSA:  3/11  implicants  -  0.27 


(EXAMPLE  5)  USING  DUECK  AND  MILLER  OM  tt2.mvl,  OUTPUT  THE  ORIGINAL 

EXPRESSION  AND  K-MAP  WITH  FINAL  RESULT. 

(command  Line  is  as  follows . ) 

ece: /work/mag 
%  mvl  -Dem  tt2.mvl 

(the  following  appears  on  thet  screen . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 
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Expression : 


radix:  4 

nvars:  2 

nterms:  11 
coeff:  2 
XI  ( 1  ,3) 

X2 ( 1 , 1 ) 
coe  f  f :  1 
XI ( 1 , 3 ) 

X2( 1 , 3) 
coeff:  2 
Xl(l, 2) 

X  2  ( 0 , 2  ) 

COe  f  f  :  1 

XI  ( 1 , 2 ) 

X2( 3, 3) 

coe  f  f :  2 
XI ( 0 , 2 ) 

X2( 0, 2 ) 
coeff:  2 
XI ( 3, 3 ) 

X2( 2,3) 
coeff:  3 
XI ( 0 , 2 ) 

X2 ( 0, 0 ) 
coeff:  1 
XI  ( 1 , 2 ) 

X2 ( 0, 2 ) 
coe  f  f :  1 
XI  (  0 , 3  ) 

X2( 1 , 2 ) 
coeff:  1 
XI ( 1 , 2) 

X2( 1 , 2 ) 
coeff:  1 
XI ( 1 , 2 ) 

X2  ( 0 , 1 ) 

Orig  map  (D6M): 

3.  3.  3.  0 

3.  3.  3.  3. 

3.  3.  3.  3. 

0  2  2  3. 

3.  3.  3.  0 

3.  3.  3.  3. 

3.  3.  3.  3. 

0  0  0  3. 

3.  3.  3.  0 

3.  3.  3.  4. 

3.  3.  3.  4. 

0  0  0  4. 

4.  4.  4.  0 

4.  4.  4.  4. 

4.  4.  4.  4. 

0  0  0  4. 

D&M:  3/11  impllcants  -  0.27 
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(EXAMPLE  6)  USING  GOLD  ALGORITHM  ON  tt2.mvl,  OUTPUT  THE  ORIGINAL  EXPRESSION 
K-MAP  AND  EACH  IMPLICANT  FOUND  WITH  FINAL  RESULT. 

(  command  line  Is  as  follows . ) 

ece: /work/mag 
%  mvl  -Geim  tt2.mvl 

(  The  following  appears  on  the  screen . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 


Expression ! 


radix;  4 
nvars;  2 
nterms;  11 
coe f ft  2 
XI  ( 1 , 3 ) 
X2(l,l) 
coeff;  1 
XI  ( 1 , 3 ) . 
X2( 1 , 3 ) 
coeff:  2 
XI ( 1 , 2 ) 
X2( 0 , 2 ) 
coe  f  f :  1 
XI ( 1 , 2 ) 
X2 ( 3 , 3 ) 
coe  f  f :  2 
XI  ( 0 , 2 ) 
X2(0, 2 ) 
coeff:  2 
XI ( 3 , 3  ) 
X2  ( 2 , 3  ) 
coe  f  f :  3 
XI ( 0/ 2) 
X2 ( 0 , 0 ) 
coe  f  f :  1 
XI  ( 1 , 2 ) 
X2 ( 0, 2 ) 
coeff:  1 
XI  ( 0 , 3  ) 
X2  ( 1 , 2 ) 
coe  f  f :  1 
XI ( 1 , 2 ) 

X2 ( 1 , 2 ) 
coe  f  f :  1 
XI ( 1 , 2 ) 
X2(0, 1 ) 


Orig  map  (D&M) : 


3. 

3. 

3. 

0 

3. 

3. 

3 

3. 

3. 

3. 

3. 

3. 

0 

2 

2 

3. 

DsM 

MIM : 

2*X1( 

Imp: 

2‘X1( 

3. 

3. 

3. 

0 

3. 

3. 

3. 

3. 

3. 

3. 

3. 

3. 

0 

0 

0 

3. 

1) *X2(  3, 

2 )  *X2(  3, 


3) 

3) 
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D&M 

HIM 

3  *  XI  ( 

Imp 

3*X'  ( 

3. 

3. 

3 

.  0 

3. 

3. 

3, 

4  . 

3. 

3. 

3, 

4  . 

0 

0 

0 

4. 

D&M 

MIM: 

3  *X1  ( 

Imp: 

3  ‘XI  ( 

4  . 

4  . 

4  . 

0 

4  . 

4  . 

4  , 

4  . 

4  . 

4  . 

4  . 

4  . 

0 

0 

0 

4  . 

Expression : 

radix:  4 

nvars:  2 

nterms:  il 
coeff:  2 
>51(1,3) 
>52(1,1) 
coe  f  f :  1 
>51(1  ,3) 
>52(1,3) 
coe  f  f :  2 
>51(1,2) 
X2(0,2) 
coeff:  1 
>51(1,2) 
>52(3,3) 
coeff:  2 
XI ( 0 , 2 ) 
X2(0, 2) 
coe  f  f :  2 
XI ( 3 , 3 ) 
X2 ( 2 , 3 ) 
coe  f  f :  3 
X1{0, 2) 
X2 ( 0, 0 ) 
coeff:  1 
XI ( 1 , 2 ) 

X2 ( 0 , 2 ) 
coe  f  f :  1 
XI ( 0 , 3  ) 

X2 ( 1 , 2 ) 
coe  f  f :  1 
XI ( 1 , 2) 

X  2  ( 1 , 2 ) 
coeff:  1 
>51(1,2) 
X2( 0, 1 ) 


Orig 

map 

(PSA): 

3. 

3. 

3.  0 

3. 

3. 

3.  3. 

3. 

3. 

3.  3. 

0 

2 

2  3. 

P&A 

MIM: 

3  *X1  ( 

Imp: 

3*X1( 

,  3 ) *X2(  3,  3) 
,  3 ) *  X2 (  1,  3) 


,  2 ) *X2(  0,  0) 
,  2 ) *X2 (  0,  2) 


1 )  *X2(  a,  o) 

2) ‘X2(  0,  2) 
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4. 

4.  4. 

0 

4. 

4.  4  . 

3. 

4  . 

4  .  4  . 

3. 

0 

2  2 

3. 

PSA 

MIM: 

3  *  XI  ( 

3, 

3)*X2( 

3,  3) 

Imp: 

3  *X1  ( 

3, 

3)*X2( 

1-  3) 

4. 

4  .  4  . 

0 

4  . 

4  .  4  . 

4  . 

4  . 

4.  4. 

4. 

0 

2  2 

4. 

PSA 

MIM: 

2  *X1  ( 

1, 

1)*::2< 

3,  3) 

Imp: 

2  *X1  ( 

1, 

2)‘X2( 

0,  3) 

4. 

4.  4  . 

0 

4. 

4  .  4  . 

4  . 

4. 

4  .  4  . 

4  . 

0 

0  0 

4. 

Gold  (  "=  )  : 

3/11 

implicants 

(Mote/  Gold(-— )  means  that  both  PsA  and  D&M  have  the  same  result  and 
the  resulting  expression  of  DsM  is  chosen  on  ties.) 


(EXAMPLE  7)  USING  GOLD  ON  tt2.mvl,  OUTPUT  STATISTICS  ON  SELECTED  ALGORITHM 

WITH  FINAL  RESULf. 

(command  line  is  as  follows . ) 

ece: /work/mag 
%  mvl  -Gs  tt2.mvl 

(The  following  appears  on  the  screen . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 

Gold(-— ).*  3/11  implicants  -  0.27 

Statistics  for  DSM: 

Num  Num  Num  Avg  Hum 

Radix  Var  Expr  Terms  Term/Fxpr 

42  1  11  11.00 

Eval  Comp  Pick 

Expr  CF  MIM 

Tot:  1254  173  4 

Avg/Expr :  1254.00  173.00  4.00 

Avg/Term:  114.00  15.73  0.36 

Statistics  for  PsA: 

Num  Num  Num  Avg  Num 

Radix  Var  Expr  Terms  Term/Expr 

42  1  11  11.00 

Eval  Pick  Gen  Next 

Expr  MIM  Bounds  Impl 
Tot:  534  4  3  30 

Avg/Expr:  534.00  4.00  3.00  30.00 

Avg/Term:  48.55  0.36  0.27  2.73 


Valid 

Impl 

27 

27.00 

2.45 


Comp  Gen 

RBC  Bounds 

23  3 

23.00  3.00 

2.09  0.27 


Next  Valid 

Impl  Impl 

27  24 

27.00  24.00 

2.45  2.18 
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(EXAMPLE  8)  USING  GOLD  ON  tt4.mvl,  VERIFY  THE  MINIMAL  SOLUTION  BY  TEST 

VECTORS  SPECIFIED  IN  THE  INPUT  FILE. 

(command  line  is  as  follows . ) 

ece: /work/mag 
%  mvl  -Gv  tt4.mvl 


(  the  following  appears  on  the  screen . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 
Gold(— —  )t  3/11  implicants  -  0.27 

Verification : 


Xl<0,3) 

X2( 0 , 3 ) 

(0,0)  -  3 

(1,0)  -  3 

(2,0)  =  3 

(3,0)  -  0 

(0,1)  -  3 

(1.1)  -  3 

(2.1)  -  3 

(3.1)  -  3 

(0,2)  -  3 

(1.2)  -  3 

(2.2)  -  3 

(3.2)  -  3 

(0,3)  -  0 

(1.3)  -  2 

(2.3)  -  2 

(3.3)  -  3 


(EXAMPLE  9)  USING  GOLD  ON  tt5 .mvl, OUTPUT  THE  ORIGINAL  EXPRESSION, IMPLICANT 
FOUND,  K-MAP,  STATISTICS  ON  THE  SELECTED  ALGORITHM , AND  VERIFI¬ 
CATION  OF  FINAL  RESULT. 

INPUT  EXPERSSION  HAS  A  DON’T  CARE  TERM  IN  IT. 

(command  line  is  as  follows . ) 

ece: /work/mag 
%  mvl  -Geimsv  tt5.mvl 

(The  following  appears  on  the  screen . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 


Expression: 

radix:  4 

nvars:  2 
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nterms:  11 
coeff:  2 
Xl(l,3) 

X2( 1, 1 ) 
coeff:  1 
XI ( 1  (  3 ) 
X2(l, 3) 
coeff:  2 
Xl( 1, 2) 
X2(0, 2) 
coeff:  4 
X 1  ( 1 , 2 ) 

X2 ( 3 , 3 ) 
coeff:  2 
XI ( 0 , 2) 

X2 ( 0 , 2 ) 
coe  f  f :  2 
XI  ( 3 , 3 ) 

X2( 2 , 3 ) 
coe  f  f :  3 
XI  ( 0 , 2 ) 

X2(0, 0) 
coeff:  1 
Xl( 1, 2) 

X2 ( 0 , 2 ) 
coeff:  1 
XI ( 0, 3 ) 

X2 ( 1 , 2) 
coeff:  1 
XI ( 1 , 2) 

X2( 1 , 2 ) 
coeff:  1 
XI ( 1 , 2) 

X2( 0, 1 ) 

Orlg  map  (DsM) : 


3. 

3. 

3 

0 

3. 

3. 

3 

.  3. 

3. 

3. 

3 

3. 

0 

4 

4 

3. 

DSM 

MIM: 

3 *X1 (  3 

Imp: 

3*XI (  C 

3. 

3. 

3. 

0 

4. 

4  . 

4. 

4  . 

4. 

4  . 

4  . 

4  . 

0 

4 

4 

3. 

DSM 

MIM: 

3  *X1 (  3 

Imp: 

3*X1(  3 

3. 

3. 

3. 

0 

4. 

4  . 

4  . 

4  . 

4. 

4  . 

4  . 

4. 

0 

4 

4 

4. 

DSM  M1M:  3*X1(  2 
Imp:  3*X1(  0 

4.  4.  4.  0 

4.  4.  4.  4. 

4.  4.  4.  4. 

0  4  4  4. 


,  3 ) *X2 (  1,  1) 
,  3 ) *X2 (  1,  2) 


/  3 ) * X2(  3,  3) 
,  3 ) *X2 (  3,  3) 


,  2) *X2(  0,  0) 
,  2) *X2(  0,  0) 
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Expression : 


radix:  4 

nvars:  2 

nterms:  11 
coe  f  f  :  2 
Xl ( 1 , 3 ) 

X  2  ( 1 , 1  > 
coe  f  f :  1 

Xl( 1, 3 ) 
X2( 1, 3) 
coeff:  2 

XI  ( 1 , 2 ) 
X2(0, 2) 

coe  f  f :  4 
Xl ( 1 , 2 ) 
X2(3,3) 
coeff:  2 
X1(0, 2) 
X2(0, 2) 
coeff:  2 
Xl ( 3 , 3 ) 
X2( 2 , 3 ) 
coeff:  3 
X1(0, 2) 
X2( 0>  0) 
coeff:  1 
X 1  ( 1 , 2 ) 
X2{ 0 , 2 ) 
coe  f  f :  1 
Xl ( 0 , 3 ) 
X2( 1 , 2 ) 
coeff:  1 
Xl(l,2) 
X2( 1 , 2 ) 
coeff:  1 
Xl ( 1 / 2 ) 
X2( 0, 1 ) 


Orig  map  (PSA): 

3.  3.  3.  0 

3.  3.  3.  3. 

3.  3.  3.  3. 

0  4  4  3. 

P6A  HIM:  3 *X1 {  2,  2)*X2(  1,  lj 
Imp.  3  *X1 (  0,  2 ) *X2 (  0,  2) 

4.  4.  4.  0 

4.  4.  4.  3. 

4.  4.  4.  3. 

0  4  4  3. 

PSA  MIM:  3 *X1 (  3,  3)*X2(  2,  2) 
Imp.  3*X1<  1,  3 ) *X2 (  1,  3) 

4.  4.  4.  0 

4.  4.  4.  4. 

4.  4.  4.  4. 

0  4  4  4. 


Gold(PsA)!  3/11  implicants  - 
Verification:  , 

Xl(0,3) 


0.27 
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X2(0, 

(0,0) 

3) 

3 

(1,0) 

- 

3 

(2,0) 

- 

3 

(3,0) 

- 

0 

(0,1) 

- 

3 

(1,1) 

- 

3 

(2,1) 

- 

3 

(3,1) 

- 

3 

(0,2) 

- 

3 

(1,2) 

3 

(2,2) 

- 

3 

(3,2) 

- 

3 

(0,3) 

- 

0 

(1,3) 

- 

3 

(2,3) 

- 

3 

(3,3) 

- 

3 

Statistics 

for  D&M 

Num  Num 
Radix  Var  Expr 

4  2  1 


Num  Avg  Num 

Terms  Term/Expr 
11  11.00 


Tots 

Avg/Expr: 
Avg/Term ! 


Eval  Comp  Pick 
Expr  CF  M1M 

1334  169  4 

1334.00  169.00  4.00 

121.27  15.36  0.36 


Statistics  for  P&A: 


Num  Num 
Radix  Var  Expr 

4  2  1 


Num  Avg  Num 

Terms  Term/Expr 
11  11.00 


Comp  Gen 

RBC  Bounds 

31  3 

31.00  3.00 

2.82  0.27 


Next 

Impl 

36 

36.00 

3.27 


Eval 

Pick 

Gen 

Next 

Valid 

Expr 

MIM 

Bounds 

Impl 

Impl 

Tot: 

811 

3 

2 

54 

52 

Avg/Expr : 

811.00 

3.00 

2.00 

54.00 

52.00 

Avg/Term: 

73.73 

0.27 

0.18 

4.91 

4.73 

(EXAMPLE  10)  USING  GOLD  ON  tt2.mvl,  IF  USER  DOESN'T  NEED  FINAL  RESULT 
SPECIFY  -q  FLAG. 

(command  line  is  as  follows . ) 

ece: /work/mag 
*  mvl  -Gq  tt2.mvl 

(  the  following  appears  on  the  screen . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 


» 

(d)  HOW  TO  GENERATE  LAYOUT  FILE  IN  MAGIC  FORMAT. 


Valid 

Impl 

33 

33 . 00 
3.00 
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SPECIFY  INPUT  ,  OUTPUT  FILE  AND  THE  SELECTED  MINIMIZATION  ALGORITHM. 
THE  FINAL  RESULT  WILL  BE  STORED  IN  THE  OUTPUT  FILE  WITH  FORMATTED 
SOLUTION. 

(command  line  Is  as  follows,  input  file:  tt2.mvl,  outfile:  tt2.out 
and  output  file  should  have  -o  flag.) 
ece:/work/mag 

%  mvl  -G  tt2.mvl  -ott2.out 

(  The  formatted  final  solution  will  be  saved  in  an  output  file  which 
will  be  used  for  layout  generation.  If  the  minimization  result  is 
worse  than  the  original  expression,  the  original  expression  will  be 
output  for  producing  the  layout.) 

(  the  following  appears  on  the  screen . ) 

MVL  VI. 0  Copyright  198B  Naval  Postgraduate  School 

Gold(— —  )i  3/11  implicants  -  0.27 


(The  formatted  output  is  as  follows.) 

I  of  input  I  of  output 
coeff 

lower  upper 
lower  upper 

99  (sentinel) 

/*  I  of  input  -  1  of  variables 

I  of  output  -  I  of  expressions  */ 


(an  example  of  the  formatted  output) 

2  1 
2 

1  2 
3  3 
3 

3  3 
1  3 
3 

0  2 
0  2 
99 

(  The  following  command  initiates  the  layout  generation 
by  using  output  file  , called  tt2.out  ) 
ece; /work/mag 
%  mvpla  tt2.out 


(The  following  appears  on  the  screen  . ) 


Input  i  tt2.out 

Output  file  i  tt2.out.mag 
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Check : 

2  inputs  1  outputs  function 


fO  -  2x0(l,2)xl(3,3) 
+  3x0( 3, 3 )xl ( 1 , 3  ) 
+  3x0(0, 2 )xl ( 0 , 2 ) 


Now,  generating  CMOS-MV-PLA  layout  in  magic  format 

fO  -  2x0 ( 1 , 2 )xl (  3 , 3 ) 

+  3x0 ( 3 , 3 ) xl ( 1 , 3 ) 

+  3x0(0, 2 )xl ( 0 , 2 ) 

Output  to  tt2 .out. mag . . . 


Done  ! 

*************  EXIT  ************* 


(  Once  layout  file  with  the  Magic  format  is  created, 
the  user  needs  to  invoke  the  Magic  graphic  package. ) 

(  command  line  is  as  follows . ) 

ece : /work/mag 
%  Magic  tt2.out 


(The  following  appears  on  the  screen  . ) 

Magic  -  Version  4.10  -  Last  updated  12/6/85  at  18:33:06 
Using  technology  "scmos". 
tt2.out:  500  rects 


(After  the  whole  layout  is  displayed  on  the  CRT,  the  user  can  create  a  file 
for  use  in  generating  a  hard  copy  of  the  layout  and  exit  by  typing  ttie 
the  following  command  to  create  the  .cif  file  for  plotting.) 

:clf  write  tt2.out 
>  tquit 


(The  following  command  line  allows  the  user  to  have  a  hard  copy  of  layout.) 
ece : /work/mag 
%  Cif plot  -r  tt2.out.cif 
(1]  20267 

ece : /work/mag 

%  Window:  -10950  12600  -2100  51000 
Scale:  1  micron  is  0.0198588  inches  (504x) 

The  plot  will  be  0.39  feet 

(1J  Done  cif plot  -I  -P  /tools/berk86/lib/patterns  -r  tt2.out. 


(the  following  is  an  example  of  a  hard  copy  of  MVL-PLA  layout.) 
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THIS  COMPLETES  THE  DESIGN  PROCEDURE  FOR  THIS  CAD  TOOL. 


2. ANALYSIS  PROCEDURES 

THERE  IS  A  RANDOM  EXPRESSION  GENERATOR  IN 
A  SET  OF  RANDOM  EXPRESSIONS  OF  A  SPECIFIC 


THIS  CAD.  THE  USER  CAN  GENERATE 
CLASS  AND  ANALYZE  THE  DIFFERENT 


MINIMIZATION  PROGRAMS 


(A)  HOW  TO  GENERATE  RANDOM  EXPRESSIONS 

(EXAMPLE  1)  THERE  IS  HELP  MANUAL  FOR  USERS  WHICH  CAN  BE  OBTAINED  AS  FOLLOWS. 

ece: /work/mag 
%  mvltest  — 

(The  following  HELP  messages  appears  on  the  screen . ) 

MVLTEST  VI. 0  Copyright  1988  Naval  Postgraduate  School 

usage:  mvltest  (-sN  -rN  -IN  -tN  -oNJ 

-  sN  Cjolc  the  random  seen  n  -imes  (default  -  0) 

-rN  -  Set  radix  to  N  (default  -  4) 

-IN  -  Set  number  of  inputs  to  N  (default  -  2) 

-tN  -  Set  number  of  terms  to  N  (default  -  1) 

-oN  -  Set  number  of  output  to  N  (default  -  1) 


(EXAMPLE  2)  THE  FOLLOWING  COMMAND  INITIATES  A  RANDOM  FUNCTION  GENERATOR 

TO  CREATE  A  SET  OF  RANDOM  EXPRESSIONS.  IT  CYCLES  THE  RANDOM  SEED  50 

TIMES  AND  GENERATE:  20  EXPRESSIONS  OF  4-VALUED  2-VARIABLE  LOGIC  EXP 

RESSION.  ALL  OF  THEM  ARE  STORED  IN  ttS.mvl  AND  HAVE  10  PRODUCT  TERMS. 

THE  NUMBER  OF  OUTPUTS  IS  THE  NUMBER  OF  FUNCTIONS  IN  THE  SAMPLE  SET. 

(command  line  is  as  follows.'...) 
ece: /work/mag 

%  mvltest  -s50  -r4  -12  -tlO  -o20  >  tt6.mvl 


(The  following  is  a  set  of  random  expressions  generated  by 
the  previous  command. ) 

4:2: 

+2*X1 ( 1 , 3 ) *X2( 3,3) 

+1*X1 (0,0) *X2( 1,2) 

+1*X1(0,3)*X2(1,2) 

+1*X1(1,2)*X2(1,2) 

+1*X1(1,2)*X2(0/1) 

+2*X1(0,0)‘X2(1,1) 

+3*X1 ( 0, 1 ) *X2( 1 , 3 ) 

+2*X1 (3,3)*X2(0,1) 

+3*X1(0,3)*X2(1,1) 

+2*X1(0,3)*X2(0,0); 

4:2: 

+1*X1(0,3)*X2(0,3) 

+2‘X1(1,1)*X2(0,3) 

+1‘X1(1,2)*X2(0,1) 

+1*X1 (1, 1)*X2<0, 2) 

+3*X1(0,0)*X2(2,2) 

+1*X1(3,3)*X2(0,3) 

+2*X1(1,2)*X2(1,1) 

+2*X1(1,1)*X2(2,2) 

+3*X1(3,3)*X2(0,3) 

+1*X1(1,1)*X2(1,1) / 

4:2: 

+2*X1(1,2)*X2(0,1)  , 

+3‘X1(1,2)»X2(0,3) 

+3*X1(2,3)*X2(0,1) 
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+2«X1(0,1)‘X2(2,3) 

+  2 *X1 ( 1 , 3 ) *X2 ( 2 , 2  ) 

+ 1  *  XI ( 0 , 2 ) *  X2( 3 , 3 ) 
+3*X1(1,1)‘X2(0,1) 

+  1 *X3 ( 1 , 2 ) *  X2 ( 0 , 0 ) 
+2*X1(0,0)‘X2(0,3) 
+3*X1(3,3)*X2<0,2); 
4:2: 

+2*X1(1,2)*X2(3,3) 

+2*X1(1,2)*X2(2,3) 

+1‘X1(0,1)*X2(1,2) 

+1*X1(1,3)*X2(3,3) 

+3*X1(0,3)*X2(1,1) 

+1*X1(3,3)*X2(2,3) 

+3*X1(3,3)*X2(1,2) 

+  1*X1 (1,2)*X2(3, 3) 
+1*X1(0,1 )‘X2(3,3) 

+  3*Xl(0,l)*X2(l,3)  ; 
4:2: 

+2-Xl(3,3)»X2<2,2) 

+3*X1(2,2)*X2(3,3) 

+1*X1(0,1)‘X2(1,3) 

+3*X1(1,3)*X2(2,3) 

+2*X1(1,1)*X2(0,2) 

+3*X1(0,2)*X2(2,3) 

+  1*X1 (2,2)*X2(1,1) 
+1*X1(2,2)*X2(1,2) 
+1*X1(2,2)*X2(2,2) 
+3*X1(0,1)*X2(1,2) ; 
4:2: 

+3*X1(1,2)*X2(1,2) 

+3*X1(2,3)*X2(1,1) 

+1*X1(2,2)*X2(0,0) 

+3*X1(1,2)*X2(1,1) 

+1*X1(3,3)*X2(1,2) 

+3*X1(1,3)*X2(0,0) 

+2‘X1(0,2)*X2(1,2> 

+  1*X1 (0,1 )*X2(1,3) 
+1*X1(0,1)‘X2(0,0) 
+3*X1(3,3)‘X2(0,0) ; 
4:2: 

+3*X1(0,3)*X2(1,3) 

+3*X1(1,1)*X2(1,3) 

+2*X1(0,3)*X2(0,2) 

+3*XJ(2,3)*X2(0,1) 

+2‘X1(0,1)*X2(3,3) 

+3*X1(0,2)*X2(3,3) 

+3*X1(1,1)*X2(0,3) 

+3*X1(1,1)*X2(1,1) 

+2*X1(0,1)*X2(1,1) 

+  3*X1(2,2)‘X2(0,3)  i 
4:2: 

+3*X1(1,2)*X2(1,2) 

+1*X1(0,3)*X2(2,3) 

+2*X1(3,3)*X2(2,3) 

+3*X1(0,3)*X2(1,1) 

+2‘X1(2,3)*X2(2,3) 

+2*X1(0,0)*X2(0,0) 

+3*X1(0,0)*X2(1,1) 

+3*X1(0,1)*X2(1,2) 

+3*X1(2,3)‘X2(0,0) 

+3*X1(2,3)*X3(0,1), 

4:2: 

+3»X1(2,3)‘X2(0,3) 

+3‘X1(0,2)‘X2(0,2) 

+3*X1(1,2)*X2(2,2) 
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+3*X3(0,0)*X2<2,3) 
+1*X1(0,3)*X2(2,2) 
+2‘X3(1,3)*X2(0,0) 
+2*X1(1, 1)*X2(1,3) 

+ 1  *  X 1 (0, 2)*X2(0, 1 ) 
+3*X1(0,1)*X2(3, 3) 
+2* XI (2, 2)* X2 (1,1) ; 
4  :  2  f 

+  3  *XJ {2,2) *X2{  0,1) 
+1*X1(0,2)*X2(2,2) 
+3*X3 (0,3)*X2(2,3) 
+2*X1(0,1)*X2(1,3) 
+2*X1(0,0)*X2(1,2) 
+1*X1(3,3)*X2(0,3) 

+  2 * XI ( 1 , 3 ) *X2 ( 2 , 3  ) 
+1‘X1 (0, 3)*X2(0,0) 
+2*X1(3,3)*X2(1,3) 
+3*X1(0,0)*X2(1,1); 
4:2: 

+  1*X) ( 0 , 3 ) *X2 ( 2 , 2 ) 

+  3  *X1 ( 1 , 3 ) *  X2 ( 0 , 3) 
+1*X1(2,2)*X2(0,2) 

+  3  *X1 ( 1 , 1 ) *X2 ( 0 , 1 ) 
+1‘X](1, 3 ) *X2 ( 1 , 1 ) 
+1*X1(3,3)*X2(2,3) 
+3‘X1(2,3)*X2(2,3) 
+2*X1(0,1)*X2(1,1) 
+1*X1(1,3)*X2(0,2) 
+1*X1(0,0)*X2(3,3)  ; 
4:2: 

+1*X1(0,1)*X2(2,3) 

+2*X1(1,3)*X2(C,2) 

+3‘X1(1,3)*X2(2,3) 

+3 *X1(1,3)‘X2(2,3) 
+1*X] <0,2)*X2(0,1) 
+3*X1(0,1)*X2(0,0) 
+2*X1(2,2)*X2(1,3) 
+3*X1(2,2)*X2(0,3) 
U*X1U,3)*X2(0,0) 
+1*X1(0,1)*X2(0,2); 
4:2: 

+1*X1(0,3)*X2(3,3) 

+2’X1(0,3)‘X2(2,2) 

+1*X1(0,3)*X2(0,3) 

+1*X1(0,1)*X2(0,2) 

+1*X1(2,3)*X2(0,3) 

+3*X1(1,3)*X2(0,2) 

+2»X1(0,1)*X2(0,2) 

+1*X1(2,2)‘X2(3,3) 

+1*X1(1,1)»X2(0,3) 

+3*X1(1,2)*X2(2,2); 

4:2: 

+2*X1(0,0)*X2(0,0) 
+1*X1(2,3)*X2(1,3) 
+3*X1(2,3)*X2(0,3) 
+1‘X1(3,3)*X2(0,3) 
+2*X1(3,3)*X2(1,2) 
+3‘X1(0,0) *X2(2, 2) 
+2‘X1(0,1)*X2(0,2) 
+2*X1 ( 2, 3 ) *X2( 1,3) 
+2*X1 ( 1 , 1 ) *X2 ( 1 , 2 ) 
+2*X1(0,1)*X2(3,3)  i 
4:2: 

+2*X1 ( 1 , 3 ) *X2( 1 , 1 ) 
+2*X1(0,2)‘X2(0,3) 
+1*X1(0,3)*X2<3,3) 
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+  2  *  X 1 (0, 3) *X2(0,0) 

+  3*X3 (0, 2) *X2(P,  3) 
+1*X3(1,3)‘X2(1,1) 
+1‘X1(0,0)*X2(1,3) 

+  3*X3 ( 2, 2) *X2(0, 0) 
+3*X1 ( 2, 3) ‘X2(l , 2) 

+  2*X1(0,3) *X2 (0 , 2 )  ; 
4:2: 

+2*X1 ( 2, 3) ‘X2( 2, 2) 
+2*X1 ( 1/ 2) *X2(1, 1) 
+3‘XI(1,3)*X2(1,2) 
+2‘X3(3,3)‘X2(1,1) 
+1*X1(1,1)*X2(3,3) 

4  3‘X3 (0, 3) ‘X2(0,  1 ) 

+  3*X3(2,3)*X2(3.,3) 
43  *X1 ( 3, 3) *X2(2, 2) 
+1‘X3(1,2)*X2(1,3) 
+2»X1( 2, 2) *X2(2, 3) ; 
4:2: 

+3*X1(1,2)*X2(0,0) 
4-3*Xl(0,l)*X2(3,3) 
+3*X1(2,3)*X2(2,2) 
+2*X1(3,3)*X2(3,3) 
4-3.  * Xl  ( 0 , 0  )  *X2 ( 1 , 3  ) 
+1*X1(1,1)*X2(2,2) 

+  1‘X3  (0,3)*X2(1,3) 

+  3  *X1  ( 0 , 3. )  *  X2  ( 1 ,  3. ) 
+3*X1(0,3) *X2(1, 2) 
4-2*Xl(3,3)*X2(2,2); 
4:2: 

+3*Xl(0,3)*X2(0/0) 
+1*X1(0,2)*X2(1,1) 
+3*X1(0,1)*X2(2,2) 
+1*X1( 2, 2) »X2( 3 , 3 ) 
+2*X1(3,3)*X2(1,1) 
+3‘X3(0,3)*X2(1,2) 
4-1  *X1  ( 3  ,  3. )  *  X2  (  2 , 3 ) 
+3»X1(1,1)*X2(0,2) 
+3*X1(1,3)«X2(0,3) 

+  1*X3.(3,3)*X2(2,3) ; 
4:2: 

4-3*Xl(0,3)*X2(0,0) 
+1*X1(1,3)*X2(0,2) 
4-3*X3(2,3)*X2(0,l) 
+2*X3(0,1)*X2(1,1) 
+1*X1(0,2)*X2(0,0) 
+1*X1(1,1)*X2(0,1) 
+2*X3(3,3)*X2(3,3) 
+2*X1(0,3)*X2(3,3) 
+3*X1(1,3)*X2(1,1) 
+2*X1(0,2)*X2(0,3) ; 
4:2: 

+3*X1(2,2)*X2(1,3) 
+1*X3(2,3)*X2(2,2) 
+1*X1(0/0)*X2(1,2) 
4-2*Xl(0,l)*X2(0,3) 
+3*X1(0,3)*X2(0,0) 
+2*X1(1,1)‘X2(1,2) 
4-3*Xl(0,3)*X2(0,2) 
+3*X1 (2,3) *X2(0,0) 
+1*X1(0,2)*X2(0,2) 
+2*X1(1,3)*X2(0,3); 
(data  file/  tt6.mvl) 
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( EXAMPLE  3)  THE  FOLLOWING  COMMAND  LINE  APPLIES  Gold  ON  tt6.mvl.  THE 

STATISTICS  AND  THE  FINAL  RESULTS  ARE  FOLLOWED. 

ece : /woi k/mag 
%  mvl  -Gs  tt6.mvl 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 


Gold( PSA)  ! 

6/10 

implicants 

- 

0.60 

Gold<-— )  : 

6/10 

implicants 

- 

0.60 

Gold ( --- )  : 

3/10 

implicants 

- 

0. 30 

Gold( - )  ; 

5/10 

implicants 

- 

0. 50 

Gold( --- )  i 

4/10 

implicants 

- 

0.40 

Gold(DSM) ! 

4/10 

implicants 

- 

0.40 

Gold( - )  ! 

3/10 

implicants 

- 

0. 30 

Gold( - )  i 

4/10 

implicants 

- 

0.40 

Gold(” — ) : 

1/10 

implicants 

- 

0.10 

Gold ( )  * 

7/10 

implicants 

- 

0.70 

Gold  ( —  )  : 

3/10 

implicants 

- 

0.30 

Gold(PSA) i 

6/10 

implicants 

- 

0.60 

Gold ( - — ) ! 

3/10 

implicants 

- 

0. 30 

Gold ( ) : 

5/10 

implicants 

- 

0.50 

Gold(  ■ - )  j 

3/10 

Implicants 

- 

0.30 

Gold(-”“ ) t 

4/10 

implicants 

- 

0.  40 

Gold(DsM) : 

4/10 

implicants 

- 

0.40 

Gold( -*■- ) t 

2/10 

implicants 

- 

0.20 

Gold ( --- ) ! 

5/10 

implicants 

- 

0.50 

Gold(--~) : 

3/10 

implicants 

- 

0. 30 

Statistics  for  D&M: 


Num  Num  Num  Avg  Num 

Radix  Var  Expr  Terms  Term/Expr 
4  2  20  200  10.00 


Eval 

Comp 

Pick 

Comp 

Gen 

Next 

Expr 

CF 

MIM 

RBC 

Bounds 

Impl 

Tot! 

31979 

3012 

101 

943 

81 

904 

Avg/Expr ! 

1598.95 

150.60 

5.05 

47.15 

4.05 

45.20 

Avg/Term: 

159.89 

15.06 

0.51 

4.71 

0.40 

4.52 

Statistics 

for  PSA: 

.  1 

Num  Num  Num  Avg  Num 

Radix  Var  Expr  Terms  Term/Expr 


Valid 
Impl 
970 
48.50 
4 . 85 


45 


4 


2 


20 


200 


10.00 


Tot: 

Avg/Expr : 
Avg/Term: 


Eval 
Expr 
16630 
831 . 90 
83.19 


Pick 
MIM 
101 
5  05 
0.  51 


Gen 

Bounds 

81 

4 . 05 
0.40 


Next 
Impl 
1032 
51 . 60 
5.16 


Valid 

Impl 

1033 

51.65 

5.17 


( EXAMPLE  4)  HERE  GENERATE  ANOTHER  SET  OF'  RANDOM  EXPRESS IONS , 10  EXPRESSIONS 
OF  4-VALUED  3-VARIABLE  LOGIC  EXPRESSION  AND  SAVE  THEM  IN  THE  DATA  FILE 
tt7.mvl. 

(command  line  is  as  follows . ) 

ece : /work/mag 

%  mvltest  -r4  -13  -t5  -olO  )  tt7.mvl 


(EXAMPLE  5)  APPLY  Gold  OH  tt7.mvl(  GET  THE  FINAL  RESULTS  AND  STATISTICS. 

(commmand  line  is  as  follows . ) 

ece: /work/mag 
%  mvl  -sG  tt7.mvl 

(the  following  appears  on  the  screen . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 


Gold ( ) 

:  4/7 

implicants 

-  0.57 

Gold(DSM) 

:  6/7 

implicants 

-  0.86 

Gold(~”= ) 

«  7/7 

implicants 

-1.00 

Gold(PSA) 

:  8/7 

implicants 

-  1.14 

Gold(DRM) 

:  7/7 

implicants 

-  1.00 

Gold(D6M) 

*  5/7 

implicants 

-  0.71 

Gold(DSM) 

:  6/7 

implicants 

-0.86 

Gold( D6M) 

:  7/7 

implicants 

-1.00 

Gold(— “) 

«  7/7 

implicants 

-1.00 

Gold  ( ) 

7/7 

implicants 

-1.00 

Statistics 

for  D&M: 

Num 

Num 

Num  Avg 

Num 

Radix  Var 

Expr 

Terms  Term/Expr 

4  3 

10 

70  7. 

00 

Eval 

Comp 

Pick 

Comp 

Gen 

Next 

Valid 

Expr 

CF 

MIM 

RDC 

Bounds 

Impl 

Impl 

Tot: 

54192 

2525 

74 

1126 

64 

97.2 

1  257 

Avg/Expr : 

5419.20 

252.50 

7.40 

112.60 

6.40 

92.20 

125.70 

Avg/Term: 

774.17 

36.07 

1.06 

16.09 

0.91 

13.17 

17 . 96 

4f> 


Statistics  for  P&A: 


Num  Num  Mum  Avg  Num 

Radix  Var  Expr  Terms  Term, Expr 
43  10  70  7.00 


Eval 

Pick 

Gen 

Next 

Val  id 

Expr 

M1M 

Bounds 

Impl 

Impl 

Tot: 

25434 

79 

69 

1269 

1333 

Avg/Expr : 

2543.40 

7.90 

6  .  90 

126.90 

133.30 

Avg/Term: 

363.34 

1.13 

0.99 

18.13 

19.04 

(EXAMPLE  6)  OUTPUT  SOURCE  EXPRESSIONS  FOR  ALL  SOLUTIONS  WHOSE  RATIO  TO  THE 

ORIGINAL  NUMBER  OF  TERMS  EXCEEDS  1.00.  ALL  WILL  BE  STORED  IN  tt7.out 

WHICH  HAS  -X  FLAG  ON  IT. 

(command  line  Is  as  follows . ) 

ece : /work/mag 

%  mvl  -1.0  -xtt7.out  -Gs  tt7.mvl 

(the  following  appears  on  the  screen . ) 

MVL  VI. 0  Copyright  1988  Naval  Postgraduate  School 


Gold( — -) : 

4/7 

implicants  - 

0.  57 

Gold(DSM) : 

6/7 

implicants  - 

0.86 

Gold(  — )  t 

7/7 

implicants  - 

1.00 

Gold (P&A) : 

0/7 

implicants  - 

1.14 

Gold(D&M) : 

7/7 

implicants  - 

1.00 

Gold(D6M) : 

5/7 

implicants  - 

0.71 

Gold(D&M) : 

6/7 

implicants  - 

0.86 

Gold(D&M) : 

7/7 

implicants  - 

1.00 

Gold( - ) : 

7/7 

implicants  - 

1.00 

Gold(--™ ) : 

7/7 

implicants  - 

1.00 

Statistics  for 

D&M: 

Num  Num  Num  Avg  Num 

Radix  Var  Expr  Terms  Term/Expr 
43  10  70  7.00 


Eval 

Comp 

Pick 

Comp 

Gen 

Next 

Valid 

Expr 

CF 

M1M 

RBC 

Bounds 

Impl 

Impl 

Tot: 

54192 

2525 

74 

1126 

64 

922 

1  257 

Avg/Expr: 

5419.20 

252.50 

7.40 

112.60 

6.40 

92.20 

125.70 

Avg/Term: 

774 . 17 

36.07 

1.06 

16.09 

0.91 

13.17 

17.96 

Statistics  for  P&A: 


Num  Num  Num  Avg  Num 
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Radix  Var  Expr  Terms  Term/Expr 
43  10  70  7.00 


Eval 

Pick 

Gen 

Next 

Val  i  d 

Expr 

M1M 

Bounds 

Impl 

Impl 

Tot: 

25434 

79 

69 

1269 

1333 

Avg/Expr : 

2543.40 

7.90 

6.90 

126.90 

133.30 

Avg/Term: 

363.34 

1.13 

0.99 

18.13 

19.04 

(EXAMPLE  7)  THE  FOLLOWING  SHOWS  ALL  THE  EXPRESSIONS  WHICH  EXCEEDS  THE 
SPECIFIED  RATIO  IN  THE  PREVIOUS  COMMAND  LINE. 


I  PSAs  8/7 

4:3: 

+2*xl(l,2)*x2(l,2)*x3(0,l) 

+2*xl(0,l)*x2(l,3)*x3(l,3) 
+I*xl(0,0)*x2(l,l)*x3(l,l) 
+3*xl(l,3)*x2(l,3)*x3(0,3) 
+I*xl(l,l)*x2(0, . ) *x3 (0 , 0 ) 
+l*xl (0, 3 )*x2(0, 3 ) *x3(0, 2  ) 

+ l*xl( 1, 1 ) *x2( 0, 3 ) *x3 ( 0, 1 ) ; 

I  DSM:  8/7 

4:3: 

+I*xl(l,l)*x2(0, 1)*X3(2, 3) 

+2‘xl(l, 3)*x2(2, 2)‘x3<0, 3) 
+3*xl(0,2)*x2(0,l)*x3(l, 1 ) 
+2*xl(0,2)*x2(0,0)*x3(0,3) 
+3‘xl(3,3)*x2(0,2)*x3(0,l) 
+3*xl(3,3)*x2(l,l)*x3(2,2) 
+2*xl(2,3)*x2(l,l)*x3(li3); 

I  PSA:  8/7 

4:3: 

+2*xl(0,3)*x2(0, 3)*x3(l,l) 

+I*xl(3,3)*x2(2,3)*x3(l,3) 
+3*xl(l,2)*x2(0,l)*x3(2,2) 
+3*xl(l,l)*x2(3,3)*x3(2,2) 
+3*xl(2,2)*x2(3,3)*x3(0,3) 
+2*xl( 0, 1) *x2( 1, 3)*x3(l, 1 ) ; 


I  PSA:  8/7 
4:3: 

+3*xl(0,2)*x2(3,3)*x3(3,3) 

+I*xl(l,l)*x2(0,3)‘x3(0,3) 

+2*xl(l,l)*x2(l/l)*x3(lr2) 

+I*xl(0,l)*x2(l,l)*x3(3,3) 

+2*xl(0,l)*x2(2,3)*x3(l,l) 

+2*xl(0,l)*x2(0,l)*x3(0,0) 

+3*xl(2,3)*x2(l,3)*x3(0,2); 


(data  file;  tt7.out) 

i 
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APPENDIX  B 

Program  Listing 


THE  FOLLOWING  ARE  THE  MODULES  USED  FOR  THE  MVL-CAD.  THE  ENTIRE  PROGRAM 
IS  ABOUT  4,000  LINES  LONG. 

HERE  IS  A  BRIEF  OVERVIEW  OF  WHAT  IS  INCLUDED. 

(1)  VERSION  BANNER,  HELP  DISPLAY  AND  GLOBAL  DATA  STRUCTURES  FOR  MVL-CAD. 

(2)  main.c 

THIS  MODULE  1)  PROCESSES  COMMAND  LINE  ARGUMENTS  2)  OPENS  FILES  AND 
3)  CONTROLS  THE  PARSING,  LOGIC  MINIMIZATION  AND  OUTPUT  FUNCTIONS. 
ALSO,  IT  INCLUDES  Gold( ) ,  WHICH  IMPLEMENTS  THE  GOLD  HEURISTIC  ( 

WHICH  CHOOSES  BETWEEN  POMPER-ARMSTRONG  AND  DUECK-MILLER, PREFERRING 
DUECK-MILLER  ON  TIES.) 

(3)  defs.c 

THIS  MODULE  INCLUDES  ALL  GLOBAL  VARIABLE  DEFINITIONS. 

(4)  parse, 1 

THIS  MODULE  IS  A  LEX  SPECIFICATION  FILE  FOR  THE  INPUT  SCANNER/PARSER 
FOR  THE  MVL-CAD.  IT  SPECIFIES  THE  GRAMMAR  OF  EXPRESSIONS  INPUT  BY  THE 
USER.  ALSO,  IT  INCLUDES  verify(),  WHICH  ALLOWS  THE  USER  TO  DETERMINE 
THE  FUNCTION  VALUE  STORED. 

(5)  dm.c 

THIS  MODULE  IMPLEMENTS  THE  DUECK-MILLER  HEURISTIC  FOR  MINIMIZING  AN 
EXPRESSION. 

THE  FOLLOWING  SUBROUTINS  ARE  INCLUDED  IN  THIS  MODULE. 

a)  Dueck_Miller( ) 

~  performs  the  Dueck  and  Miller  heuristic  on  the  Input  expression 
and  returns  the  number  of  implicants  used  to  cover  it. 

b)  _c f ( E , X ) 

-  computes  the  clustering  factors  at  X,  where  X  is  a  vector  of 
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coordinates  of  a  non-zero  minterm  which  has  the  smallest  value. 

-  returns  the  clustering  factor,  an  integer. 

c)  *mim(E) 

-  finds  the  most  isolated  minterm  determined  by  the  minterm  with 
smallest  CF  and  returns  a  vector  of  integers  representing  the 
coordinate  of  it  or  NULL  if  no  more  minterms 

-  the  function  value  at  that  location  is  also  returned  as  the  last 
integer  in  the  vector. 

d)  valid_implicant( I ) 

-  checks  implicant  validity  and  returns  1  if  the  proposed  impli- 
cant  is  valid;  else  returns  0.{An  implicant  is  valid  if  it  is 
indeed  an  implicant  of  the  function.) 

e )  compu  te_r be ( I ) 

-  accepts  implicant  I  and  calculates  the  relative  break  count  of  it. 

f)  *pick_implicant(^' 

-  finds  the  most  isolated  minterm  according  to  the  smallest  CF. 

-  generates  all  implicants  that  covers  this  minterm. 

-  computes  the  rbc  for  each  and  chooses  the  implicant  with  the 
smallest  rbc  (best  implicant). 

-  returns  the  best  implicant. 

( 6 )  pa . c 

THIS  MODULE  IMPLEMENTS  THE  POMPER-ARMSTRONG  HEURISTIC  FOR  MINIMIZING 

AN  EXPRESSION. 

THE  FOLLOWING  SUBROUTINES  ARE  INCLUDED  IN  THIS  MODULE. 

a)  Pomper_Armstrong( ) 

-  performs  the  Pomper-Armstrong  heuristic  on  the  input  expression 
and  returns  the  number  of  implicants  used  to  cover  it. 

b)  *mim(E) 

-  finds  the  most  isolated  minterm  in  the  input  expression  based 
on  the  random  weight  generated  by  random  number  generator. 

-  returns  a  vector  of  integers  representing  the  coordinates  of 
the  most  isolated  minterm  or  NULL  if  no  more  minterms. 
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-  the  value  of  the  function  corresponding  to  the  most  isolated 
minterm  is  also  returned  as  the  last  integer  in  the  vector. 

c)  random( ) 

-  returns  a  random  number  using  a  linear  congruential  method. 

d)  valid_implicant( I ) 

-  checks  implicant  validity  and  if  valid,  returns  1  and  the  num_ 
zeroed,  the  number  of  minterms  driven  to  zero  by  this  implicant 
when  subtracted,  else  returns  0. 

e)  *pick_implicant(X) 

-  finds  the  most  isolated  minterm  randomly  and  then  produces  all 
implicants  that  cover  this  minterm. 

-  computes  num_covered,  the  number  of  minterms  covered  by  this 
implicant. 

-  picks  the  best  implicant  which  drives  the  maximum  number  of 
minterms  to  zero.  If  there  are  more  than  one  such  implicant 
the  largest  implicant  is  chosen. 

-  returns  the  best  implicant. 

( 7 )  common . c 

THIS  MODULE  INCLUDES  HEURISTIC  SUFPORT  FUNCTIONS  COMMON  TO  DUECK-MILLER 
AND  POMPER-ARMSTRONG . 

THE  FOLLOWING  SUBROUTINES  ARE  INCLUDED  IN  THIS  MODULE. 

a)  *_eval(E,X) 

-  evaluates  the  expression  at  (xl , x2, x3 xn ) 

b)  *next_coord(T) 

-  computes  the  next  possible  coordinate  for  term. 

-  returns  an  integer  vector (X)  containing  the  coordinates. 

c)  *gen_bounds(X) 

-  generates  the  permissible  bounds  around  location  X  in  the  work¬ 
ing  expression  for  use  in  picking  the  next  implicant  that  covers 
X. 

-  returns  a  bounds  array. 

d)  *next_implicant(B) 
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-  chooses  X  as  the  first  implicant. 

-  On  each  call,  returns  the  next  implicant  with  a  bound  array (B) 
from  *gen_bounds( )  that  specifies  the  implicant. 

e)  subtract_implicant( I ) 

-  adds  the  best  implicant(I)  to  the  working  expression  as  a  neg¬ 
ative  term  (negated  coef ficient) . 

-  Also  adds  the  best  implicant  to  the  final  expression. 


(8)  alloc. c 

THIS  MODULE  INCLUDES  ALL  MEMORY  ALLOCATION  FUNCTIONS.,  SUCH  AS  alloc_ 
term( ) ,  alloc_bound( ) , and  dealloc( ) ,etc.  IT  IS  USED  TO  ALLOCATE  SPACE 
FOR  A  TERM  ARRAY,  BOUND  ENTRIES  AND  DEALLOCATE  OR  DUPULICATE  THE  EXP¬ 
RESSION. 

(9)  mvltest.c 

THIS  MODULE  IS  A  TEST  SOURCE  FILE  GENERATOR.  IT  GENERATES  ANY  TYPE  OF 
MVL  EXPRESSIONS  ACCORDING  TO  THE  SPECIFICATION  GIVEN  BY  THE  USER. 

(10)  Makefile 

THIS  MODULE  KEEPS  TRACK  OF  WHICH  SOURCE  FILES  ARE  DEPENDENT  ON  OTHERS 
AND  PROVIDES  FOR  AN  EFFICIENT  RECOMPILATION  OF  A  PROGRAM. 
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/ *  mvl 


-  A  multi-valued  logic  expression  compiler/analy zer/optimiier 

LCDR  John  M.  Yurchak 
LT  Hoon  Seop  Lee 
Dr.  Jon  T.  Butler 

Copyright  1988  U.S.  Naval  Postgraduate  School,  Monterey,  California 
All  rights  reserved  */ 


f include  "defs.h" 
f include  <ctype.h> 


/*  Version  banner  and  help  display  */ 


static  char 

•version  -  "\nMVL  VI. 0  Copyright  1980  Haval  Postgraduate  School\n\n", 
•usage  -  "usage:  mvl  (-elmrsvqDFG)  (-0.0)  |-orlleJ  (-xrile)  (file 
•helpH  -  ( 


-e 

-1 

-m 

-V 

-S 

-q 

-r 

-0.0 


"  -D 

«  -p 

"  -G 

"  -xrile 

. 

"  -orlle 

HULL 


-  Print  the  original  expression,  as  parsed", 

-  Print  each  implicant  found", 

-  Print  the  Karnaugh  map  of  the  original  and  each", 
successive  expression", 

HAM1IHG:  Impractical  for  nvar  >  3", 

-  verify  the  minimal  solution", 

-  Print  statistics  on  selected  heuristics", 

-  Quiet,  don't  print  final  results", 

-  Permit  negative  coefficients  in  terms", 

-  Build  source  expressions  for  all  solutions", 

"hose  ratio  to  the  original  number  of  terms", 
exceeds  this  number  and  output  them  to  a  file", 

(default  is  \"x.mvl\"(", 

-  Execute  the  Buecfc  £  Miller  heuristic  (default)", 

-  Execute  the  Pomper  *  Armstrong  heuristic", 

-  Execute  the  Gold  heuristic", 

-  Output  source  expressions  built  using  -0.0  to  \"rlle\"", 
Instead  of  \"x.mvl\""# 

-  Ouput  formatted  solutions  to  \"rile\"". 
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/*  Global  data  structures 


/*  Logic  expressions : 

E_orig 

-  holds  the  original  input  expression  as  parsed 


B_work 

-  a  copy  a  Eorig 

-  lmplicants  are  subtracted  from  this  expression  as  terms 
during  the  coures  of  optimisation 


Efinal ( I 

-  the  result  expression  (starts  out  empty) 

-  each  term  is  one  impllcant  found  during  optlmlratlon 

-  each  heuristic  has  its  own  E_final  (for  comparison) 


Expression 

E_orig  -  (  0,0,0,  HULL  |, 
e_»ork  -  (  0,0,0, HULL  ), 
E_final(Z)  -  ( 

(  0,0,0,  HULL  |, 

(  0,0,0,  HULL  | 

)  ; 


lnt  HEUR; 


int  FIHAL? 


/*  Current  heuristic 

*  IfEUR  indexes  into  E  final  (J 

*  depending  upon  the  currently 

*  active  heuristic 
*/ 

/•  Index  of  the  selected  final 

*  expression 
*/ 


char 

expression |MAX_BUr+l), 
i  f  _name  |MAX__rATH+l  |  , 
xf_name (MAX_FATH  + 1 )  -  "x.mvl\0", 
«tmp_name  -  "mvlXXXXXX", 
ofname (MAX  PATH+l ) ; 


/*  Parser  working  buffer  */ 

/*  Input  file  name  */ 

/*  Expression  output  file  name  */ 

/*  Output  file  name  (for  result)  */ 


int  token_type; 
int  completed; 


jmp_buf  eof_context; 


/•  Scanner  token  code  */ 

/*  rlag: 

*  1  if  a  complete  sentence  has  been 

*  recognised,  0  otherwise 
*/ 

/*  Saved  context  for  lex  end  of  file  */ 


/*  user  options  */ 


int  r_flag  -  0, 
e_flag  “  0, 
iflag  -  0, 
m_flag  -  0, 
a_flag  -  0, 
v_flag  «•  0, 
D_flag  »  1, 
Pflag  «  0, 
G_flag  -  0. 
q_flag  -  0; 

int  nvar, 
radix; 


/*  Humber  of  variables  in  the  current  expression  */ 
/*  Radix  of  the  current  expression  * / 


r  * 

O'l 


int  tot_expr  -  0,  /*  Total  number  of  expressions  optimized  */ 

tot_terms  —  0;  /*  Total  number  of  terms  in  all  input  expressions 

MVL_stats 

*STAT,  /*  Pointer  to  the  current  statistics  record  */ 

OM^statf  / *  Stats  record  for  the  D8M  heuristic  */ 

FA_stat;  /*  Stats  record  for  the  P4A  heuristic  */ 

float  TO_ratio  -  0.0; 

rn.E  *xfile? 


-  yywrap ()  is  called  by  lex  on  EOF . 

-  Thus,  we  execute  a  non-local  goto  back  to  the  setjmp  call  which 
saved  the  context,  at  which  point  we  open  the  next  input  file. 

-  Hotice  that  yywrapO  never  actually  returns. 

*/ 

yywrapO 

l 

long jmp (eof_context, 1;  ? 


main (argc, argv) 
char  *argv( | 


: f unct ion  ; 

-  Process  command  line  arguments 

-  Open  files 

-  Control  the  parsing,  optimization  and  output  functions 


register  ap; 
char  *p; 

fprlntf (stderr, "%s",  version) ; 

/*  Do  the  options  */ 

for  fap-1;  ap  <  argc;  ap++)  | 

p  -  argv |ap) ; 
if  <*p41  —  '-')  | 

/*  loop  through  an  option  string  */ 
while  <*p)  I 

switch  (*p44)  ( 

case  ' D'  : 

D_f lag44; 
break; 
case  ' P’ ! 

P_f lag44; 
if  (Dflag— 1) 

P  flag--; 
break; 
case  ' G' : 

G_f lag44; 

Df lag44 ; 

P_f lag44; 
break; 
case  ' v'  : 

v  flag44; 
break; 
case  ' x' : 

if  <‘r)  I 

strcpyfxf  name,p); 

/*  terminate  the  while  loop  */ 
*P  _  -\0'  ; 


else  | 

fprint f (stderr, "mvl :  filename  expected  after  -x\n") ; 
fprintf(stderr,"%s", usage); 
exit  (1) ; 

I 

break; 
case  ' o' : 

if  <*P)  ( 

strcpyfof  name,p); 

/*  terminate  the  while  loop  */ 

*p  -  'NO'; 

I 

else  t 

fprint f  (stderr, "mvl :  filename  expected  after  -o\n"); 
fprint f (stderr, "Is",  usage) ; 
exit  (1)  ; 

I 

break; 
case  ’ g' : 

q_f lag44; 
break; 


5' 


case  ' s'  : 

s_f lag+t; 
break,* 
case  '  r'  : 

r_£lag++; 
break; 
case  ' m' : 

m_£lagt +; 
break; 
case  ' e'  : 

e_f lag++; 
break; 
case  '  i'  : 

i_f  lagf-f ; 
break ; 
case  '  , '  : 
case  ' 0' ; 
case  '  1'  : 
case  ' 2'  ; 
case  ' 3'  : 
case  ' 4'  : 
case  ' 5'  : 
case  '  6'  : 
case  '  T  : 
case  ' 0'  ; 
case  ' 9'  : 

P— ; 

sscanf (p,  "%f  ",  fcFO^ratio) ; 

*p  -  '\0'  ; 
break; 
case  '  : 

lprint£(stderr,"%s", usage) ; 
for  (ap-0;  help(ap);  ap+  +  ) 

fprintffstderr, "%s\n", help  lap]); 

exit  (0)  ; 
default ; 

fprintf (stderr, "mvl!  unknown  switch  -  %c\n" , p ( -1 ) ) ; 
f print  £ (stderr, "%s", usage) ; 
exit  (1) ; 


I 

) 

I 

/*  Erase  the  old  expression  file  iif  there  is  one)  */ 
unlink (xf_name) ; 

/*  If  the  user  specified  an  output  file  ...  */ 
if  <of_name 10] )  ( 

/*  erase  the  existing  output  file  and  open  the  new  one  (append)  */ 
uni Ink (of_name) ; 
irk  temp  (tmp_name)  : 

if  ( (yyout”£open (tmp_name, "a") )  MULL)  ] 

fprintf (stderr, "mvl :  can't  open  temp  file  -  %s\n", tmp_name) ; 
exit (1) ; 

I 

I 

/*  Init  the  scanner  */ 
token_type  -  TOK_EOL; 

/*  For  each  input  file  */ 
for  (ap-1;  ap  <  argc;  apt+)  ( 
p  -  argv(ap) ; 

/*  skip  the  option  arguments  */ 
if  <*p  —  '-') 
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cont Inuf ; 

strcpy  (if_name, p)  : 

if  (  (yy  irv»£open  ( i  f  _name,  "r  "  )  )  NULL)  ( 

fprintf (stderr, "mvl :  can't  open  %s\n", if_name) ; 
exit  (1)  ; 

I 

inlt_9tats (SPMstat)  ; 
init  9tata (SFA_stat) ; 


/*  Set  the  EOr  context  */ 
if  (set jmp (eof_context)  —  0)  ( 


/*  Hhlle  9t i 1 1  expte99ion9  to  paree  In  this  file  ...  */ 
while  (expr()  >  0)  ( 

nvar  —  E  orlg.nvar; 
radix  —  E_orig . radix; 
tot_terms  4™  E_orig . nte tm; 
tot  _expr44; 


I 


if  (D_flag)  (  /*  do  the  Dueck  S  Miller  */ 

Dueck_Mlller  O ; 
riHAL  -  D_M; 

I 

if  (Pflag)  )  /*  do  the  Pomper  t  Armstrong  */ 

Fomper_Armst rong  ()  ; 
if  (I0_flag) 

riNAL  -  PA; 

I 

if  (Gflaq)  I  /*  do  the  Gold  (which  chooses  between  above) 
Gold!)  ; 

if  (of_name(0))  /*  output  the  final  expression  */ 
output  expr (HEUR) ; 
riNAL  -  HEUR; 

) 

else  1 

if  (of_name|0))  t 

if  (D_£lag)  /*  output  the  final  expression  */ 
outputexpr (D_M) ; 
else  if  (P_f lag) 

output_expr (P_A) ; 

I 


I 


verify))  ; 

dealloc_expr (SE_orig) ; 


I 

else  if  ( I  completed) 

syntax_error (ERR_UNEXP_EOf) ; 


*/ 


/*  Print  the  statistics  */ 
if  (s  flag)  ( 
if  (Dflsg) 

print_DM_stats  ()  ; 
if  (P_f  lag) 

print_FA  stats  (); 

» 

f  close  (yyin)  ; 

) 

if  (of  name(OJ)  ( 


/*  Build  the  output  file  for  the  layout  generator  by  appending 

*  the  list  of  final  expressions  In  the  temporary  file  to 

*  a  header  consisting  of  radix,  inputs  and  outputs. 

* 

*  final  format  is: 
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radix  inputs  (nvar)  outputs (tot^expr) 
coef  f 

lower  upper 
lower  upper 

99  (sentinel) 
coef  f 

lower  upper 
lower  upper 


*  99  (sentinel) 

* 

*  etc 

V 

/*  Close  the  temporary  file,  reopen  for  reading  and  open 

*  the  final  output. 

V 

f close (yyouh ) i 

if  (  (yyin  «■  f open  ( tmp_name  ,  " r  ** )  )  *■*  NULL)  ( 

f print f (stderr , "mvl :  can't  open  temp  file  -  %s\n" , tmp_name) ? 
exit  (1) ? 

) 

if  (  (yyout  —  fopen  (of_name,  "w")  )  —  MULT  }  ( 

fprintf (stderr, "mvl:  can't  open  output  file  -  %s\n" , of_name) ; 
exit  (1)  ; 

) 

/*  Build  and  output  the  header  */ 

sprintf (expression, "%d  %d  %d\n", radix, nvar , tot_expr) ; 
f puts (expression, yyout)  ; 

/*  Append  the  temp  file,  then  delete  it  */ 
while  (fgets (expression, MAX_BUF, yyin)  I-  MULL) 
fputs (expression, yyout ) ; 
f close (yyout ) ; 
f close (yyin) ? 
unlink (tmp_name) ; 

I 

) 
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Gold  0 

/*  - - 

: function : 

-  Implement  the  Gold  heuristic  (which  chooses  between  and  D4M, 

preferring  D6H  on  ties)  . 

: global s ! 

E_final() 

E  orlg 
HEUR 

! called_by ! 
main  () 

- - -  */ 


char  *p; 

if  (E_f inal  ID_MJ  . nterm  <  E_f inal IP_A) . nterm)  ( 

HEUR  -  P_M? 
p  -  "D&M"; 

) 

else  if  (E_f i nal |D_M) . nterm  >  E_f inal |P_A) . nterm)  ( 

HEUR  -  FA; 
p  -  "P4A"; 

) 

else  f 

HEUR  -  D_M; 

p  »  **  "  • 

i 

if  (Iqfl aq)  I 

prlntf("  Gold(Va):  %4d/%-4d  implicants  -  %4 . 2£%%\n\n", 
p,  E_f  inal  (HEt’R|  .  nte tm,  E  or  ig  .  nterm, 

(float) E_f inal (MEUR) . nterm/ (float) E_orig . nterm)  ; 

t 

) 


print_map { ) 

/*  - 

: f  unct ion : 

-  Feint  the  Karnaugh  map  of  E_work  in  lta  present  state 
tglobals : 

E_work 
:called_hy : 

Dueck_Mi Her  ( ) 

Potnper_Armstrong  O 
:  calls : 

eval  0 


( 

register  i, j; 
int  X(MAX_VAR+1 )  ; 
int  *  V; 

for  (1-0;  i  <  nvar;  1  +  +  )  X [ i J  -  0; 
for  (i-0;  i  <  nvar;)  I 

V  “  ^eval  (re  work,X); 

printf  ("»3d%c",  VlEVAh] ,  V|I1LV]  '  )  ; 

x  ( i  I  +  +  ; 

for  (;i  <  nvar;)  ( 

if  (X[i)  >-  radix)  f 
X(i)  -  0; 
if  (i  <  2) 

printf ("\n") ; 

111; 

x  ( i  J  +  + ; 

I 

else  ( 

i  -  0; 
break; 


print 

implicant (s,  X,  I ) 

char 

*  s; 

int 

*X; 

Term 

*1? 

/  — 

function; 

-  Print  the  Host  Isolated  Hinterm  X  and  the  implicant  selected 
to  cover  it  I . 

*  cal led_by : 

Dueck_Millet  {) 

Pomper_Armst rong  () 


register  i; 

printf  ("  %s  HIM;  %2d", s, X { nvar J ) ; 
for  (i-'O;  i  <  nvar;  144) 

printf ( " *X%d (%2d,  %2d)  ”,±  +  1,  X(1),X(1M  ; 
printf ("\n") ; 

printf  ("  Imp:  %2d" , I->coef f ) ; 

for  (i*=0;  i  <  nv?4r;  144) 

printf ("*X%d  (%2d,  %2d) ", i+l, I->B[i) . lower, I->B(i) .upper) ; 
printf  r\n\n")  ; 

I 

print_t erms  (E) 

Expression  *E; 

/* - 

! function : 

-  Print  the  expression  pointer  to  by  E 
: calledby ; 

DueckMi ller  <) 

Poir>per_Armst  rong  () 


int  i, j; 

printf ( "\nExpression : \n\n" ) ; 
printf ("  radix:  %2d\n", E->radlx) ; 
printf i"  nvars:  %2d\n", E->nvar ) ; 
printf <"  nterms:  %2d\n" , E->nterm) ; 
for  (1--0?  i  <  E->nterm;  i44)  [ 

printf("  coeff:  %d\n" ,  E->T ( i J . coef f ) ; 
for  ( j^O;  j  <  E->nv«r;  j  +  4) 

printf  ("  X%d (%d, %d) \n", 

341, E->T ( i ) . B ( 3 1 • lower,  E->T[i) .B(3J  .upper)  ; 

) 

printf  <"\n") ; 

) 
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output_expr (heur) 
int  heur; 

/*  - 

:  function : 

-  Append  the  final  expression  to  the  temporary  file. 

-  heur  indexes  into  E  final. 

-  If  the  final  expression  is  not  at  least  as  good  as  the  original 
input  expression,  output  the  original. 

-  Output  format  as  follows: 

coef  f 

lower  upper 
lower  upper 


(sentinel) 


99 

:  globals : 
yyout 
stdout 
E_orig 
E~final 
:  called_by : 
main ( ) 


register  i,  j; 

Expression  *E; 

if  (yyout  1—  stdout)  ( 

if  (E_orig . nterm  <  E_f inal ( heur J . nt erm) 

E  -  fiE_orig; 

else 

E  -  SE_f inal [ heur J ; 
for  (i-0;  i  <  E->nterm;  i  +  +  )  ( 

f print f  (yyout, "%d\n",E->T(i| . coef f )  , 
for  (j”0;  j  <  E->nvar;  j44) 
fprintf  (yyout , "%d  %d\n", 

E->T ( i ) . B l j ) . lower, 

E->T(i) . B ( j | .upper) ; 


fprintf  (yyout, "99\n") ; 
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print ^source (E,g,n) 

Expression  *E; 
char  *s; 

/*  - 

: function : 

-  Print  the  exp  1 pointA';*  t*'  by  **.  **o  ►h*  file  xf^name  in 
the  form  recogni reable  by  the  parser. 

:globals : 

xf  name 
: called_by : 

Dueck_Miller  () 

Pomper_Armst rong () 


register  ir  j; 

FILE  *f  ? 

if  ((f  -  fopen (xf_name, "a") )  —  HULL)  ( 

f print f (stderr, "mvl :  can't  open  %s\n" ,  xf__name)  ; 
return; 

) 

f pr  int  f  (  f ,  "f  %s:  %d/%d\n" ,  s,  nf  E->nterm)  ; 

fprintf  (f , "%d : %d : ", E-> radix, E->nvar ) ; 

for  (i«0;  i  <  E->nterm;  it  +  )  { 

fprintf  (f ,  "\n  +%d", E->T |i J  .coef  f )  ; 

for  (j~0;  j  <  t ->nvar;  jit)  ( 
fprintf <f, " * x%d ( %d,  %d)  ", 

j+1 >  E->T ( i J . B ( j 1 .lower, E->Tti) .B[ j) .upper) ; 

I 

) 

fprintf (f, ";\n") ; 

f close (f ) ? 

) 


init_stat 3 (S) 

MVL_st at  3  *S; 

/* - - - 

! function : 

-  Initialize  statistics  structure  pointer  to  by  S 
: called_by : 


tot_expr  “  0; 
tot_terms  -  0; 

S->cal 1 s_eval  -  0; 

S->c«lls_cf  *0; 

S->eal ls_mim  —  0? 
S->call9_gen_bound5  -  0; 
S->calls_pick_implicant  ”  0; 
S->cal  ls_next__impl i  cant  -  0; 
9->calls_valid_implicant  "  0; 
3->calls_compute_rbc  —  0; 

I 


print  PA  stats  () 

/* - 

: function : 

-  Print  the  stats  on  the  P£A  heuristic 
! globals : 

tot_terms 
tot_expr 
radix 
nvar 
PA_stat 
:called_by : 
main  () 


printf  ("Statistics  for  P&A:\n\n"); 

printf  ("  Hum  Hum  Hum  Avg  Num\n"); 

print f ("Radix  Var  Expr  Terms  Term/Expr\n") ; 
print  f  (" %4d  %4d  %6d  %6d  %8.2f\n\n", 

radix, nvar, tot_expr, tot_terms, (float) tot_terms/ (float) tot_expr ) ; 
printf("  Eval  Fick  Gen  Next  ValidNn")? 

prlntf("  Expr  MIM  Bounds  Impl  Impl\n"); 

print  f  ( "Tot !  %11  Id  %01d  %81d  IBld  %81d\n", 

PA_stat .call s_eval , FA_stat . callsmim, 

P A_st at .  cal  1  s_gen  bounds,  FA  st. at .  calls_next_implicant , 

PA  st at . calls_val id_implicant ) ; 
printf  ("Avg/Expr:  %11.2f  %0.2f  18. 2f  %8.2f  %8.2f\n", 

(float) PA_stat . calls_eval/ (float ) tot_expr, 

(float) PA_stat .calls_mim/ (float) tot_expr, 

(float) FA_stat . calls_gen_bounds/ (float) tot_expr, 

(float) PA_stat . call s_next_impl leant / (float) tot_expr, 

(float)PA  stat. calls  valid  implicant/ (float ) tot_expr) ? 
printf ("Avg/Term:  %U.2f  *0.2?  %8.2f  %8.2f  *8.2£\n\n", 

(float) PA  stat . calls_eval/ (float) tot_terms, 

(float) PA_stat ,calls_mim/ (float) tot_terms, 

(float) FA_stat . cal ls_gen_bounds/ (float) tot^terms, 

(float) PA_stat . calls_next  impl leant/ (float) tot_terms, 

(float) PA_atat . cal ls_valid_impl leant/ (float) tot_terms) ; 
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pr int_DM_stat 3  () 

/*  - 

: function : 

-  Print  the  state  on  the  D6M  heuristic 
: globals : 

tot_terms 
tot  expr 
radix 
nvar 
OM_stat 
: cailed_by : 
main  0 


printf  ( "Statistics  for  D£M:\n\n")/ 

printf  ("  Num  Num  Hum  Avg  Num\n"); 

print f ("Radix  Var  Expr  Terms  Term/Expr \n" ) ; 
print  f  ( "%4d  %4d  %6d  %6d  %8.2f\n\n"f 

radix,  nvar,  tot__expr,  tot_terms,  (float)  tot  terms/  (float)  tot_expr)  ; 
printf  ("  Eval  Comp  Pick  Comp  G**n  Next 

printf  ("  Expr  CE  MIM  P.BC  Bounds  Impl 

print  f  ( "Tot :  %llld  Iflld  181d  181d  IBld  181d  181d\n", 

DMstat .  call  s_eval,  DMstat .  calls__cf ,  DM__stat .  calls_mim, 

DM_stat . calls_compute_rbc, 

D11_st at .  cal  ls_gen_  bounds,  DM_st at .  calls_next_impl leant, 

DMstat . cal ls_valid_ impl leant ) ; 

printf ("Avg/Expr :  %11.2f  18. 2f  18. 2f  10. 2f  18. 2f  %8.2f  %8.2f\n", 

(float)  DM_stat  .calls_eval/  (float)  tot_expr, 

(float) DM  st at . eal ls_cf / (float) tot_expr, 

(float) DM  st at . cal ls_mim/ (float) tot_expr, 

(float) DMstat .calls_compute_rbc/ (float) tot_expr, 

(float) DM_stat . calls_gen_bounds/ (float) tot^expr, 

(float) DMstat . call s_next_impl leant / (float) tot_expr, 

(float) DMstat . calls  valid_implicant/ (float ) tot  expr)? 
printf ("Avq/Term:  %11.2f  18. 2f  18. 2f  10.2f  18. 2f  %8.2f  18. 2f\n\n", 

(float) DM_stat . cal ls_eval/ (float) tot_terms, 

(float) DMstat . callage f / (float) tot^terms, 

(float) DMstat .callsmim/ (float) tot  terms, 

(float) DM  stat . ca 1 ls_compute_rbc/ (float) tot_terms, 

(float) DM_stat .call s_gen_bounds/ (float) tot^terms, 

(float) DMstat .call s_next_impl leant / (float) tot _terms, 

(float)  DM_stat .  calls_valid_implicant /  (float)  tot__terms)  ? 


fatal  (s) 
char  *s; 

l 

fprintf (stderr, "%s\n", s) ; 
exit (1) ; 

) 


V  a  1 1  d  \  n 
Impl \n" 
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/*  defs.h 


-  External  definitions  for  mvl 


*/ 


•include  <stdio.h> 
linclude  <setjmp.h> 

I define  MAXBUr  4096 
•define  MAXPATH  64 
I define  HAX_VAR  32 
Idefine  MAX_INT  32765 

Idefine  YES  1 
•define  NO  0 

Idefine  EVAL  0 
Idefine  HLV  1 

•define  DM  0 

Idefine  P_A  1 

/*  parser  definitions  and  data  structures 


•define 

TOK 

EOF 

0 

•define 

tok" 

"number 

2 

Idefine 

tok" 

"VAR 

3 

Idefine 

tok" 

"colon 

4 

•define 

TUK[ 

STAR 

5 

Idefine 

tok" 

"lfaren 

6 

Idefine 

tok’ 

COMMA 

7 

Idefine 

TOK 

"rfaren 

8 

Idefine 

tok’ 

signed 

9 

Idefine 

tok" 

"dontcare 

:  10 

Idefine 

tok; 

SEMI 

11 

Idefine 

TOK 

"plus 

12 

Idefine 

TOK 

[minus 

13 

Idefine 

H 

O 

EOL 

14 

Idefine 

TOK 

3YNERR 

-1 

Idefine 

EBP. 

RAP I X_ PNG 

-2 

•define 

EPR 

"nvar  png 

-3 

Idefine 

ERR 

"coErr  png 

-4 

Idefine 

ERR 

VARINPEX 

RTG 

Idefine 

ERR 

PUP  VAR 

-6 

Idefine 

ERR 

LOWER_RNG 

-7 

Idefine 

ERR 

UFFER_RNG 

-8 

Idefine 

ERR 

UFPERLOWER 

-9 

Idefine 

ERR 

unexp  Eor 

-10 

extern  char  expression(MAX_BUF+ll ; 
extern  int  token_type; 

extern  FILE  *yyin, *yyout ; 

/*  MVL  definitions  and  data  structures  */ 

typedef  struct  exprstruct  Expression; 
typedef  struct  term_struct  Term; 
typedef  struct  bound  struct  Bound; 
typedef  struct  term_struct  Implicant; 

struct  expr_struct  | 
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int  radix, 
nvar, 
nterm; 

Term  *T; 

In¬ 
struct  term_struct  ( 
int  coeff; 

Bound  *B; 

In¬ 


struct  bound_struct  ( 
int  lower, 
upper; 

In- 


extern  Expression 
E_or 1 g, 

E_work , 

E_final|2J ; 

extern  int  !!EUR,riNAL; 

extern  int 
nvar, 
radix; 

Term  *alloc_term|) ; 

Bound  *alloc_bound () ; 

extern  char 

ifname (MAX_P ATH4 1 J , 
xf  name  (I1AXJ ATI14 1  ]  , 
of_name(MAX_PArH41  J  ; 

extern  rflag, 
e_f lag, 
i_£ lag, 
x_£ lag, 
mf lag, 

G_f  lag, 
v_f lag, 
q__f  lag; 

extern  TILE  *xfile; 

extern  int  completed; 

extern 

int  tot  expr, 
tot_terms; 

struct  mvl^stats  ( 
long 

call s_eval, 
calls  cf, 
cal ls_mim, 
ca lls_gen_bounds, 
ca 1 ls_next_impl leant , 
call s_pick_impl leant, 
cal ls_va 11 d_impl leant , 
call  s_compu  t  e__r  be; 

I; 


typedef  struct  mvl__stats  MVL_stats 
extern  MVL  stats 
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*  STAT , 

DHstat, 

FA_9t at , 

G_9tat ; 

extern  float  FO_ratio 

lnt  *mim < ) : 

lnt  *next_coord () ; 

Bound  *gen^bounds 0 ; 

Term  *next_impl leant () ; 
Term  *pick^irnpllcant ( ) i 
int  *eval  ()",  *_eval  ( ) ; 
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/*  parse. 1 


-  lex  specification  file  for  the  input  scanner /parser  for  the  MVL 
parser /compiler 

expression  grammars 

espr  !• 

radix  9:r  num_var  ' term_list 

term_list  :*• 

term  '  !  9 

l  term  '4'  term__list 


term  i- 

coeff  9  *9  var  list 


vsr__llst  :■ 

variable 

|  variable  var_list 


variable  : — 

var__id  '  ['  limit  limit  ')' 

verivicatlon  teat  grammar: 

veri f 1 cet ion  : ” 

variable  9  ? 9 

|  variable  verification 

constraints ! 

radix  >—  2 
num_var  1 

1  <-  coeff  <-  radix  (radix  —  don't  care) 

0  O  limit  <-  radix-1  and  lower__limit  <-  upper_llmit 
"xl"  <-  var_id  <-  "x (num_var ) " 

data  structures 

expression 

|  radix  ) 

(  nvar  J 
(  nterm  ) 

|  term  *  1 

-  An  expression  structure  contains  the  radix  of  the  logic 
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t  coeff  ) 

(  bounds  *J 


-  A  term  structure  contains  a  coefficient  and  a  pointer  to  a 
list  of  bounds  (length  equal  to  the  number  of  variables 
allowed  in  the  expression) 

bounds 

(  lower  ) (  upper  ) 

(  lower  ) (  upper  J 


-  A  bounds  list  is  an  array  of  pairs  (lower  and  upper  bounds) 
which  may  range  over  the  radix  of  the  expression,  one  pair 
for  each  possible  variable  in  the  term 

V 

•include  "defs.h" 

%% 

"#" I A\n) *\n  ( 

) 

i\t  j*  i 

» 

10-9J+  ( 

strcat (expression, yytext ) ; 
return (TOK  NUMBER); 


(+-M0-91+  ( 

strcat (expression,  yytext)  ; 
return (TOK  SIGNED); 


•?"  I 

strcat (expression, yytext)  ; 
return (TOK  DONTCARE) ; 


IxXl | 1-9J |0-9) *  ( 

strcat (expression,  yytext) ; 
return (TOK_VAR) ; 


s  •  «  | 

strcat (expression, yytext)  ; 
return (TOK_COLON) ; 

I 

«*"  j 

st  rest (expression, yytext) ; 
return  (TOK  STAR); 


■<"  ( 

strcat (expression, yytext) ; 
return <TOK_LPAREN) ; 

I 

"  +  "  ( 

strcat  (expression, yytext) ; 
return  (TOK_FL(JS)  ; 

I 

"  _ «  ) 

strcat (expression, yytext) ; 
return (TOK_MINUS) ; 

I 

( 

strcat (expression, yytext) ; 
return (TOK_COMMA) ; 

I 

")"  ( 

strcat (expression, yytext) ; 
return (TOK_RFAREN) ; 

I 

.  < 

return (NO) ; 

) 

1 

strcat (expression, yytext) ; 
return <TOK_SEMI) ; 

) 

t  \t\f\n\r)  ( 

expression(O)  «*  '  \0’ ; 

) 

%% 

/*  parser  NT  functions  for  expressions  */ 

expr  0 

I 

int  rad,nvar; 

inlt_expr ( ) ; 

expression 1 0)  —  '  \0' ; 

if  ( ! yytext ( V ]  ||  match(";’)) 

next_token ( ) ; 
completed  «■  0; 
if  (token_type  «"•  TOK_EOF) 
return (0) ; 

if  ( !  radlx_val  ()  )  ( 

syntax_error  (TOK_STHERR(  ; 
return (NO) ; 

) 

sscanf (yytext, "*d", trad) ; 
if  (rad  <  2)  ( 

syntax  error (ERR_RADIX_RNG)  ; 
return (NO) ; 

) 

E_orig. radix  —  rad; 

next  token  (); 
if  ( frnatch  ("!*))  ( 
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syntax  er ror  (TOK__SYNERR)  ; 
return (NO) ; 

I 

next  token  ( )  ? 
if  (lnum_var ())  1 

syntax^error (TOK_SYNERR) ; 
return  (NO) ; 

i 

sscanf  (yytext,  "%d"/Snvar)  ; 
if  (nvar  <  1)  | 

syntax  error (ERR_NVAR_RNG) ; 
return (HO) ; 

) 

E_orig.nvar  —  nvar; 

next  token ( ) ; 
if  (Tmat  ch  (":"))  I 

syntaxerror (TOK_SYNERR) ; 
return  (NO) ; 

I 

next  token  ( ) ; 

if  (Tterm_list  ()  )  ( 

syntaxerro  -  <TOK_SYNERR)  ; 
return (NO) ; 

I 

if  I  (match  {";  ") )  ( 

syntax  error (TOK_SYNERR) ; 
return (NO) ; 

) 

completed  -  1; 
return (YES) ; 

I 

radix_val () 

1 

if  (tokentype  --  TOK_NUMBER) 
return  (YES) 
return (HO) ; 

) 

num_var  () 

( 

if  (token  type  ~  T0K_NUI1BER) 
return (YES) ; 
return (NO) ; 

) 

term_list  () 

( 

register  i; 
if  ( I  term  ()  ) 

return  (NO)  ; 

if  (match ( "  +  " )  |  match ("-") )  ( 

next_token ( ) ; 
while  (term_list  () )  ; 

) 

else 

while  (term  l 1st  ( ) ) ; 
for  (1-0;  i  <  E_orig.nvar;  lit)  ( 

if  (E_or ig . T ( E_or ig . nterm-1 | ,  B  (  i  ) .lower 
E_or ig . T ( E_or ig . nt erm-1 ) ,B(i| .lower 

I 
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return ( YES)  ; 


term  () 

( 

if  ("(token  tyre  —  TOKJltlMBER)  II  (token_type  —  T0K_S1GHED)  I  I 
jjcanf (yytext, "Id", tcoef f ) ; 
if  ( I rf lag)  ( 

if  ((coeff  <  1)  !l  (coeff  >”  Eorig. radix) >  I 
syntaxercor  <ERR_COErr_RHG)  ; 
return (NO) ; 

I 

I 

else  if  (token_type  —  TOK_DONTCARE)  ( 
coeff  -  E_orig . radix; 

I 

else  ( 

return (HO) ; 

I 

£  crlg  .  ntermM  ; 

E_orig.T  -  slloc_term(E_orig. T, coeff, E_orig . nterm) ; 

next  token  f ) ; 

if  U'match  <*•*•*)  )  ( 

syntax  error  (TOK_S  YHEF.R)  * 
return^KO) ; 

I 

next  token ( ) ? 

if  (Tvar  llstt)) 

return (HO) ; 

return (YES) ; 


var_list  U 

( 

if  ( ! vst iab.1 e  () ) 
return (HO) ; 
next_t oken  ( ) ; 
if  (match ("•"> >  f 
next  token  0 ; 
while  (var_l i at ( ) ) ; 

I 

return (YES) ; 


variable  ( ) 

( 

int  vsr index, loner, upper ; 

if  ( !  var  id  ()  ) 
ret  urn (MO) ; 

sscanf fr yytext (1), w%d", fcvarindex) ; 

if  ((vtr index  <  1)  |(  (vsrindex  >  E_or ig . nvsr )  )  ( 

syntax  error (ERR_VARINDEX_RHG) ; 
return (HO) ; 

) 

var index — ; 

if  (E  orlg. T(E_orlg . nt«rm-! ) - B I var index) . lower  I-  -1)  I 
syntax  error (ERR_DUP— VAR) ; 
return(NO) ; 

I 

next  t  ^V.en  ( )  ; 

if  (imatchC  (")  )  ( 


syntax  error  (T0K_9YNERR) ? 
return  (NO) ; 

) 

next  token ( ) ; 

if  (TlinitO)  { 

syntax__er ror  (TOK_SYNERR)  ; 
return (NO) ; 

I 

sscanf  (yytext, "%d", Slower) ; 

if  ((lower  <  0)  ||  (lower  >  E_orig . radix-1 ) )  ( 

syntax  error  (ERR__LOWER_RNG)  ; 
return (NO) : 

) 

next  token  {) ; 
if  (TmntcM",  ")  )  ( 

syntax__error  (TOK_SYNERR)  ; 
return (NO) ; 

I 

next  token  ()  ; 

if  (TllmitO)  ( 

.  syntax  error (TOK_SYHERR) ; 

return (NO) ; 

I 

sscanf (yytext, "%d", Supper) ; 

if  ((upper  <  0)  ||  (upper  >  E  or ig . radix-1 ) )  ( 

syntax_error (ERR_UFPER_RNG) J 
return (NO) ; 

I 

if  (upper  <  lower)  | 

syntax  error  (ERR_UFPER__LOWER)  ; 
return  (NO) ; 

) 

next  token  ( ) ; 
if  ( (match  (")"))  ( 

synt ax_er ror (TOK_SYNERR) ; 
return (NO) ; 

I 

E  orlg.T(E_orig.nterm-3 ) . B ( varindex) . lower  —  lower 
E_orig.T(E_orig.nterm-lJ .B(varindex) .upper  -  upper 
return (YES) ? 

) 

var_id  () 

I 

if  (tokentype  --  TOK_VAR) 
ret urn (YES) ; 
return (NO)  ; 

I 

limit  () 

( 

if  (tokentype  —  TOK_NUMBER) 
return (YES) ; 
return (NO) ; 

\ 

/*  verification  NT  functions  */ 

verify  U 

( 

regi ster  ( : 

Bound  * ( MAX_VAR+ 1 | ; 

Term  T; 


int  f ir st , f i r st_t ime  —  1; 
int  ‘value, *X; 

for  (;;)  ( 

if  (match  (";")) 

next__token  ( )  ; 

T.B  -  B; 

for  (1=0;  i  <  nvar;  i4  +  ) 

B ( i I . lower  -  -1 ; 
if  ( I vervariable (B) > 
ret  urn; 

completed  -  0; 

while  (ver _var iable <B) ) ; 

if  ( (match  C;  ")  ) 

syntax_error (TOK_SYNERR) ; 
completed  —  1; 
if  (Ivflag) 
continue; 

first  -  1; 

If  (first  _t ime )  ( 

print f ("Verification: \n\n* ) ; 
first_tlme  —  0; 

I 

for  (i=0;  i  <  nvar;  i 4  4  > 

print  f ("  X%d (%d, *d) \n", ill, B|i) . lower, B |i)  . upper) ; 
pr int f ( " \n  ; ? 

while  ( (X  “  next^coord (ST, first) )  I—  NULL)  ( 
first  —  0; 

value  ~  _eval (*E_final (riNAL) ,X) ; 
printfC  <"); 
for  (1=0;  i  <  nvar;  14+)  ( 

print f  ("%d", X ( i  1 )  ; 
if  (ill  <  nvar)  printfC,"); 

I 

printfC)  -  %3d\n", value |EVAL))  ; 

) 

prlntf ("\n") ; 

) 

I 

ver_variable (B) 

Bound  *B; 

( 

int  varindex, lower, upper; 
if  ( !  var_id ( ) ) 
return (No); 

sscanf (Syytext (1) , "%d", Svarinlex) ; 
if  ((varindex  <  1)  ((  (varindex  >  nvar))  ( 

syntaxerror (ERR_VARINDEX_RNG) : 
returnluo) ; 

I 

varindex — ; 

if  (B(varindex) . lower  I-  -1)  ( 

syntaxerror (ERR_DUP_VAR) ; 
returnTNO) ; 


next  token!); 
if  (Tmatch ("("))  ( 

syntax  error (T0K_S1NERR) ; 
return (No) ; 

I 

next  token  ()  ; 

if  (Tlimit  M )  ( 

syntax_error  (T0K__S1NERR) ; 
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return (NO) ; 


I 

sscanf  (yytext,  l,%d',/  Slower)  ; 

If  ((lower  <  0)  II  (lower  >  radix-1))  I 
sy  ntax^er  ror  (ERR_LOWER__  RNG)  ; 
return  (NO)  ; 

) 

next  token  ()  ; 

if  (imatch  (’,"))  I 

syntaxerror (TOK_SfNERR) ; 
return (NO) ; 

I 

next  token  ()  ; 

it  (TllmitO)  I 

syntax  er ror (TOK_3YNERR) ; 
return (NO) ; 

) 

sscanf (yytext, "%d". Supper) ; 

if  ((upper  <  0)  II  (upper  >  radix-1))  ( 
syntax_error (ERR^UFPER_RNG) ; 
return (NO) ; 

I 

if  (upper  <  lower)  ( 

syntaxerror (ERR_UFFER_LOWER) ; 
return (NO) ; 

I 

next  token ( ) ; 

if  (Imatch  (")") )  I 

syntax  error (TOK_SYNERR) ; 
return (NO) ; 

I 

B I varindex) . lower  -  lower; 

B I varindexl . upper  -  upper; 

next_token  I) ; 

return (KES)  ; 


/*  scanner  utility  functions - - -  *i 

syntax_error (code) 

I 


static 

char  *etrors(J  • 

i 

/* 

TOK_SYNERR 

*/ 

i»  •* 

"syntax  error", 

/* 

ERR  PADIX^FNG 

♦/ 

"radix  out  of  range  (>“  2) ", 

/* 

EPF  NVAR_RNG 

♦/ 

"1  of  vat  "ut  of  range  (>-  l)"r 

/* 

EPP  COEFTRNG 

*/ 

"coeff  out  of  range  (1  <■  coeff 

<-  radJ x-1 ) ", 

f 

E  R  R_  V  A  R I  MD  EX_RNG 

*/ 

"var  Index  out  of  range  (1  to  # 

of  var) ", 

/* 

ERR_DUP_VAR 

*/ 

"duplicate  variable  id", 

l* 

EPP  LOWERING 

*/ 

"lower  bound  out  of  range  <0  to 

radl x-1 ) " , 

/* 

ERR  UrPERRNG 

♦/ 

"upper  bound  out  of  range  (0  to 

radlx-1 ) ", 

/* 

EPP  UrrER_LOWER 

*/ 

"upper  bound  must  be  >■*  lower  bound". 

/*• 

ERP_UHEXP_EOr 

*/ 

"unexpected  end  of  file". 

); 

NULL 

code  “ 

abs (code) ; 

pr  int f  ("%s\nls\n", errors (code) , expression); 
exit  (1) ; 
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next  token!) 


( 

tokentype  -  yylex(); 
if  <token_type  — «  TOK_SYHEPR)  ( 
syntax_error (TOK_SYNERR) ; 

I 

r»  urn (token_type) ; 

) 

match (s) 
char  *si 
I 

if  ( 1 strcmp (s, yytext) ) 
return (YES) ; 

else 

return (MO) ; 
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•s 


/ *  dm.  c 


-  This  module  implements  the  Dueck  t  Miller  heuristic  for  minimizing 
an  expression. 


•include  "defs.h" 

Dueck_Miller () 

/* - 1 - 

! function : 

-  Perform  the  Dueck  &  Miller  heuristic  on  the  input  expression 
j  algoritlim: 

Start  with  a  working  copy  Rwotk  of  the  original  function  E_orig; 
Initialize  a  filial  function  E  final; 

While  (there  are  still  minterms  to  pick)  ( 

Pick  a  minterm  X  from  Rework; 

Pick  the  best  implicant  1  for  X; 

Subtract  1  from  Rework; 

Add  I  to  E__f inal; 

) 

jglobals i 
R_or ig 
e_f  J  ag 
mf  lag 
qflag 
G_f lag 
TOrat io 
! side_ef f ects ! 

STAT 
IIRUR 
Ework 
E_£inal l J 
Scalled^by  j 
main  () 
icalls: 

dealloc_expr  () 
dup^expr  () 
pr int_terms () 
pr  int_map  () 
mim  () 

pick_lmplicant () 
subt ract_implicant () 
print_source ( ) 


int  num  “  0; 

int  *X; 

Term  *1; 
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float  ratio; 

if  ( E_£ 1 na 1 ( D_M]  .T  !-  NULL) 

dealloc_expr <SE_f inal (D_M) ) ; 

STAT  -  SDMstat; 

HEUR  -  DM; 

dup  expr  (*E_»*ork,  *E_orig)  ; 

Eflnal (HE0R1 .nterm  -  0; 

E_f  Inal  II1EUR)  .radix  -  E_orig  .  radix; 

E~f  Inal  ItlEUP]  .nvar  «  E_orlg.nvar; 

E_final(HEUR] .T  -  NULL; 

if  (e  flag) 

print_terma (SE_orig)  ; 
if  (m^f  lag)  ( 

prlntff"  Orig  map  (DSM):\n"); 
print  map  <)  ; 

) 

for  ( ; ; )  ( 

if  ((X  -  mim(SE_work) )  —  NULL) 
break; 

I  —  pick  impllcant (X) ; 
ntim  imp 1+4; 

subtract_implicant (I) ; 
if  (iflag) 

print  implicant ("DSM",X,  1)  ; 
if  (m_f lag) 

print_map ( ) ; 

I 

ratio  -  ( (float) numimpl/ (float) E_orlg. nterm) ; 

If  (foratio  ! “  0.0)  ( 

if  (ratio  >  ro_ratlol  ( 

pcint_source (SB_orig, "DSM", num_impl) ; 

I 

1 

if  ( I g  flag  SS  !G_flag)  I 

prlntff"  D(mT  %4d/%-4d  lmpllcants  -  M.2f*%\n\n" 
num_impl, E_orig . nterm, ratio) ; 

) 

dealloc_expr  (SE_*»ork)  ; 
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Static  int  cf(E,x) 

Expression  *E; 
int  x; 

/*  - 

: function : 

-  Determine  the  clustering  factor  for  (xt , x2, x3, x4 . . . . xn) 

-  Puilds  a  vector  of  coordinates  for  _cf() 

-  Local  to  dm.c 
:calls : 

_cf  () 

: returns : 

-  An  integer  clustering  factor 


int  *X; 

X  «  (int  *)  fix  ? 
return (_c£ (E. X)) ; 


static  int  cf(E,X) 

Expression  *E? 
int  *X; 

/* - - - 

: function : 

-  Compute  the  clustering  factor  at  X,  where  X  is  a  vector  of 
coordinates . 

-  Local  to  dm.c 
: globals : 

radix 

nvar 

: sideef f ects : 

STAT 

r called  by : 

mim  0 
: calls : 

_eval  () 
vcopy  () 

: returns : 

-  An  integer  clustering  factor 


int  nterm  -  F.->nterm; 
register  i , j, k; 
int  value ( 2 ) r 
expanded, 
va 1 1 ( 2 ) , 
val2 ( 2 ) ? 
int  cf  -  0, 
dea  —  0r 
ea  ■*  0; 

STAT->call2_cf *4 ; 
vcopy (value, _eva 1 (E, X) ) ; 

if  (  (value(RVAL)  <  1)  | |  (value(EVALj  >-  radix)) 

return (MAX_IHT) ; 

/*  for  each  variable  (direction) ...  V 
for  (i-0?  i  <  nvar;  i+4  )  ( 

expanded  —  0; 

/*  If  not  on  a  left  hand  edge,  move  left  */ 
if  <X  1 1 1 )  | 

X|  I|  —  ; 

vcopy (vail , _eval (E,  X)  )  ; 
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if  (vail  (EVAt.)  >-  value  (E VAL))  I 
expanded!!; 
eat!; 

( 

x(ijt+; 

I 

/*  if  we  didn't  start  on  a  right  hand  edge,  move  right  */ 
if  (X I i J  <  (radix-1))  ( 

XIU!!; 

vcopy(val2,  eval(E,XJ); 
if  (val2(EVAt.)  >-  value (EVAL) )  I 
expanded!! ; 
ea!!; 

I 

xii|  —  ; 

I 

if  (expanded) 
deal!; 

I 

/*  compute  the  clustering  factor  */ 

cf  —  (dea  *  (radix-1) )  1  ea; 
return (cf ) ; 


82 


Static  int  *mim(E) 

Expression  *E; 

/* - 

: function : 

-  Find  the  Most  Isolated  Minterm  in  the  expression  pointed  to 
by  E,  and  return  its  coordinates  as  a  vector. 

-  Local  to  dm.c 
tglobals : 

nvar 

:side  effects: 

STAT 

: called^by : 

Dueck_Miller  () 

: calls : 

next  coord ( ) 

_eval  () 
vcopy ( ) 

_cf  () 

:  returns : 

-  A  vector  of  integers  representing  the  coordinate  of  the  most 
isolated  minterm,  or  NULL  if  no  more  minterms. 

-  The  value  at  that  location  ia  also  returned  as  the  last  Integer 
in  the  vector. 


int  cur_val  -  E->radix, 
cur_Cr  -  MAX_INT, 
cf , 

value  12} , 
term; 

int  *X, *next_coord  () , 1 • 

static  int  save_coord [MAX_VAR41 } ; 

STAT->calls_mim4+; 

for  (term-0;  term  <  E->nterm;  termH)  ( 
i  -  l; 

while  ( (X-next_coord (& (E->T ( term) ) , i) )  I-  NULL)  ( 
vcopy  (value, _eval (E,  X) ) ; 
if  (value l EVAL J  >  0)  ( 

if  (valuefEVALJ  <  cur^val)  ( 
cur_val  -  value (EVAL); 

for  (i-0;  i  <  E->nvar;  i  +  +  )  aave_coord ( i )  -  X(i); 
curler  -  _cf  <E,  X) ; 

) 

else  If  (value(EVAL)  cur__val)  ( 
cf  -  _cf (Ef X) ; 
if  (cf  <  curCr)  ( 
cur_Cf  -  cf; 

for  (i«>0;  i  <  E->nvar;  i44)  save_coord  ( i )  -  X(i); 

) 

) 

) 

i  -  0; 

) 

) 

if  (curCr  —  MAX_INT) 
return  (NULL) ? 

save_coord(E‘>nvar|  cur_val; 
return (save_coord)  ; 

) 
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static  int  valid_implicant (I) 

Term  *1; 

/*  - 

j  function: 

-  Decide  upon  the  validity  of  implicant  I 

-  Local  to  dm.c 
: globals : 

E_wo  r  k 
B_or ig 

j  aide_ef fects : 

STAT 

: called_by : 

pick__irnpl  leant  <) 

: calls . 

next_cocrd ( ) 

_eval  () 

vcopy 0 
: returns : 

1  if  a  valid  implicant 
0  if  not 


int  *X; 

int  init  -  1; 

int  value  -  I->coeff; 

int  Vo | 2  J  r  Vw ( 2 )  ; 

STAT->calls_val id_impl leant  4  4  ; 
while  (  (X  -  next  coord  (I,  init )  )  !  -  NULL)  ( 

init  -  0; 

vcopy ( V  w  r _ e v a 1 (&E_work,  X) ) ; 

vcopy  (Vo, _eval (SE_orig,  X)  )  ; 

if  (((Vw(EVAL)  <  value)  &&  IVh(HLVJ)  &&  (Vo|EVAl)  <  (radix-1))) 
return (0) ; 

) 

r  e  t  u  r  n  ( 1 )  ; 


) 


static  int  compute_rbc { I ) 

Term  *  I ; 

/* - 

! funct ion : 

-  Compute  the  RBC  for  the  given  implicant 

-  Local  to  dm.c 
: globa 1 s : 

radix 

nvar 

:  side__ef  f  ects  : 

STAT 

j  called_by ! 

pick_implicant () 

2  cal Is: 

next_coord  0 
_eval  () 
vcopy  () 
i  returns : 

-  an  integer  RBC 


int  *X; 

int  I__value  —  I->coeff; 
register  i; 
int  val ue ( 2  J , 

neighbor_value  ( 2 ) , 

good, 

bad, 

rbc  -  0, 
init  —  1; 


STAT->calls_comput e_rbc+4; 

/+  for  each  coordinate  in  the  implicant  ...  */ 
while  ((X  -  next_coord (I, init ) )  f-  NULL)  ( 
init  -  0; 

vcopy  (value, _eval (6E_work, X) ) ; 
if  (value (EVAL)  radix) 
continue; 

/*  for  each  direction  ...  */ 
for  (i~0;  i  <  nvar;  i+4)  ( 

good  —  0 ; 
bad  -  0; 


if  (value(EVAL)  —  I_value) 
good  —  1 ; 


/*  if  there  is  a  left  neighbor,  examine  it  */ 
if  ((X(i)  I-  0)  S(  (X l i J  -«*  I->B[iJ  .lower/ )  ( 

X[iJ—  ; 

vcopy  (neighbor__value, __eval  (6E_work,  X)  )  ; 

X[il++; 

if  (neigbbor_value (EVAL)  f—  0)  ( 

if  (neighbor_value (EVAL)  —  (value (EVAL)  -  I 
good  »  1; 

if  (neighbor_value (EVAL)  value  (EVAL)) 
bad  —  1 ; 

) 

) 


/*  if  there  is  a  right  neighbor,  examine  it  V 
if  ((X(iJ  I-  (radix-1))  65  (X [ i )  —  I->B(i) .upper) ) 
X(i)+4; 

vcopy (neighbor_value,_eval (6E_work, X) ) ; 

X  ( i )  — ; 

if  (neighbor_value (EVAL)  ! ■  0)  f 

if  (neighbo revalue ( EVAL)  »  (value (EVALJ  -  I 


value)  ) 


( 


value) ) 


good  ~  1 ; 

if  (ne ighbor_value l EVAL)  — ™  value (EVALJ) 
bad  -  1; 

I 

} 

/*  update  the  rbc  */ 
tbc  -  (rbc  -  good)  +  bad; 

) 

) 

return (rbc) ; 
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static  Term  *pick_implicant (X) 
int  *X; 

/•  - - - 

! function : 

-  rick  the  beat  implicant  for  minterm  X 

: global a : 

radix 

: aide_ef  f ect a : 

STAT 

: cal led_by : 

Dueck_Miller  () 

:  calls : 

ini t_iroplicant () 
gen_bounds  0 
next  lmplicant () 

_eval  () 
vcopy ( ) 
compue^rbc ( ) 
copy_impl leant () 
va 1 id^impiicant () 

:  returns : 

-  A  pointer  to  a  term  representing  the  beat  Implicant. 


int  cur_rbc  —  MAX^INT, 
rbc  —  0; 

Implicant  *1; 

static  Bound  Inbound [MAX_VAR+ 1 J ; 

static  Term  I_best/ 

Bound  *B; 
int  V ( 2 J ; 

STAT->calls_pick_implicant+4; 

I_best.B  -  I_bound; 
init^implicant (X)  ; 

B  -  gen_bounds (X)  ; 

while  (<I  -  next_tmplicant (F) )  t-  NULL)  ( 
vcopy (V,  eval (sc  orig,X)); 
if  (V  ( EVAL J  •=-  (radix-1 )  )  ( 

for  (I->coef£  -  XfnvarJ;  I->coeff  <  radix;  (I->coef f ) ++)  ( 

if  (valid_lmpl leant (IM  f 
rbc  -  compute_rbc (I )  ; 
if  (rbc  <  cur^rbc)  ( 
cur_rbc  -  rbc; 
copy_impl leant  (fil__best, I) ; 

1 

) 


I 

else  { 

I->eoeff  "  XfnvarJ; 
if  (val4 d_implicant (I)  )  ( 

rbc  -  compute_rbc (I) ; 
if  (rbc  <  cur^rbc)  ( 
cur_rbc  —  rbc; 
copy^implicant (SI_best,  1}  i 

) 

\ 

) 

) 

return (fi I_best )  ; 

I 
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/*  pa.c 


-  This  module  implements  the  Fomper-Armst tong  heuristic  for  minimizing 
an  expression. 


•  include  "defs-h'" 

Fomper_Armst rong  O 

/*  — - - 

j  function: 

-  Perform  the  Pomper  &  Armstrong  heuristic  on  the  input  expression 
:globals : 

FO_rat io 
E_or lg 
e_f lag 
m_f lag 
q_f lag 
i_flag 
G_f lag 

: sideef f ects : 

STAT 
IIEUR 
Ework 
EfinallJ 
jcalledby : 

main  0 
scalls : 

dealloc_expr  () 
dup  expr  () 
pr int_tetms  {) 
print _map  () 
mim  () 

pick_impllcant  0 
aubtract^implicant  0 
print_source  U 


int  num_^impl  —  0; 
int  *X; 

Term  M; 
float  ratio; 

if  (E_f  inal  IP  AJ  .T  1-  MUM,) 

dealloc_expr  (*E_f  inal  (P_AJ ) ; 


STAT  -  SFAstat; 

HBUR  -  PA; 

dup^expr ( £ Rework , 4E_orig) ; 
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E_£ inal  [ IIEUR |  . nte rm  —  0; 

E_£lnal ( HEUR ]. radix  “  E_orig . radix; 

E_f 1 nal ( HEUR J , nvar  —  E_orig.nvar; 

E_£inal(HEUR)  .T  -  NULL; 

i £  <e_f lag) 

prlnt_terma (*E_orig) ; 
if  (m_f lag)  ( 

prirtf ( "  Orig  map  (PSA):\n"); 
print_map ( ) ; 

) 

for  (;;)  | 

If  MX  -  mim(«E_work)  )  --  HULL) 
break; 

I  —  pick_implicant (X) ; 
num  impl++ ; 

9ubtract_implicant  (I) ; 
if  (l_f  lag) 

print_impl leant ("PtA",  X,  1)  ; 
if  (m_£lag) 

print_map  ()  ; 

I 

ratio  -  ( <float ( num_impl/ (float) E_orig. nterm) ; 
if  (roratio  I-  0.0)  ( 

if  (ratio  >  ro_ratlo)  | 

prlnt_aource (SE_orig, "Pt A",  num_impl) ; 

1 

I 

if  <1q  flag  ««  IG  flag)  ( 

print f  ("  PSA;  »4d/»-4d  implicanta  -  %4.2f%%\n\n 
num_lmpl,  E_orig . nterm,  ratio)  ; 

) 

dealloc_expr (SE_wcrk) ; 


static  int  *mim(E) 

Expression  *E; 

/*  - 

: function! 

-  Return  a  vector  corresponding  to  the  coordinates  of  the  current 
Most  Isolated  Minterm 

•globals ! 
radix 
nvar 

: side_ef  f ect s : 

STAT 

:called_by : 

romper__Armst  rong  () 

:calls: 

next  coord O 

_eval<) 

vcopy  () 
random  () 

S  returns : 

-  A  vector  of  integers  representing  the  current  Most  Isolated 
Minterm,  or  NULL  if  no  more  minterms. 

-  The  last  integer  in  the  vector  is  the  value  of  the  function  at 
that  coordinate. 


int  weight  -  MAX_1HT; 
int  tmp; 

int  *X, 1 , *next_coord 0  ; 

Int  first; 

static  int  save_coord (MAX_VAR+1 )  -  (  -1  )? 
int  value(2); 
int  term; 

STAT->calls_mim+4 ; 

for  (term  —  0;  term  <  E->nterm;  termf+)  ( 
first  -  1; 

while  ( (X-next_coord (G (E-"T ( term) J r first) )  I-  NULL)  ( 
vcopy  (value, _e vs 1 (E,  X)  )  ; 

if  (  (value(EVAL)  >  0)  &&  (value(EVAL)  I-  radix))  ( 

tmp  —  random ( J ; 
if  (tmp  <  weight)  ( 

for  (1-0;  i  <  nvar;  i++) 
save_coord ( i J  -  X | £  J ; 
save_coord (nvar )  -  value(EVAL)? 
weight  —  tmp; 

I 

I 

first  -  0; 

I 

) 

if  (weight  1-  MAXINT) 
return (save^coord) ; 
return (NULL)  ; 

) 
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static  int  random {) 

/*  - 

: function : 

-  Return  a  number  between  0  and  32766  using  a  linear  congruential 
method 
:  cal  led__by : 
mim  () 


static  long  a— 100001; 


a  -  (a* 125)  %  2796203; 

return {  (int)  (  (  < float ) a/27  96203) *32767) ) ; 
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static  int  valid_implicant (I) 

Term  *1; 

/*  - 

: f  unct ion : 

-  Decide  upon  the  validity  of  implicant  1 
: globa Is : 
radix 

:side_ef fects: 

STAT 

:called  by: 

plck_lmplicant  () 

: cal Is : 

next_eoord ( ) 
eval  0 
vcopy  {) 

: returns : 

num_reroed  if  valid 
0  if  not 

-  v 

( 

int  ‘X; 

int  init  —  1; 

int  num_reroed  —  0: 

int  value  —  I->coeff; 

Int  Vo{2| ,V«(21 : 

STAT->calls_va 1 id_lmpl leant  4  1; 
while  ((X  -  next_coord (1,  init) )  I-  NULL)  ( 
init  *•  0; 

vcopy (Vw,_eval (SE_work,X) )  ; 
vcopy (Vo, _eval (SE_orlg,  X) )  ; 

if  (((VvtEVAL)  <  value)  CS  !Vw(HLVJ)  ll  (Vo(EVAL)  <  (radix-1))) 
return  (0) ; 

if  (Vw ( eval 1  <-  value) 
num_zeroed++; 

I 

return (num_ieroed) ; 

I 


92 


static  Term  *pick_implicant  (X) 
lnt  *X; 

/•  - 

•  function  t 

-  Pick  the  best  impllcant  for  minterm  X 
jglobals : 

nvar 

radix 

: slde_ef f ects : 

STAT 

j  called^by : 

Pompe ^Armstrong  () 

•calls : 

init_impllcant () 
gen_bounds  0 
next_implicant  () 

_eval  () 
vcopy  0 

copy_lmpllcant  () 

•  returns : 

-  A  pointer  to  a  Term  representing  the  best  impllcant. 


lnt  cur  num_zeroed  ”  -1, 
cur_num_covered  —  -1; 
lnt  num  covered, 
num zeroed; 

Impllcant  *1; 

static  Bound  Inbound (MAX_VAR4l J ? 

static  Term  I_best; 

Bound  *B; 
lnt  V  ( 2  |  ? 
register  1; 

STAT->calls_pick_impl leant  44? 

I_best.B  -  I_bound; 
init_implicant (X)  ? 

B  *=  genbounds (X) ? 

while  ((I  -  next_implicant  (P) )  I-  NULL)  ( 
vcopy (V,  eval (SE_or ig, X) ) ; 
if  (V ( EVAL  J  —  (radix-1 )  )  ( 

for  (I->coeff  -  X(nvar);  I->coeff  <  radix?  (I->coef £) 44)  ( 

if  (num_zeroed  “  valid_implicant (I) )  ( 

num__covered  ~  1? 
for  (i-0?  1  <  nvar;  144) 

nur  -overed  *«  ( I->B ( 1 ]. upper  -  1->B ( i ]. lower  4  1) ? 
if  ( (num_reroed  >  cur_num_reroed)  | |  ( 

(num_zeroed  —  cur_num_zeroed)  £4 
(num__covered  •»  cur_num_covered) 

))  i 

cur_num_reroed  «  num_zeroed; 
cur  num_covered  -  num_covered; 
copy_impl leant (Sibbest,  I)  ? 

) 

I 

I 

I 

else  ( 

I->coeff  •  X l nvar]; 

if  (num_*eroed  —  valid_implicant (I) )  ( 
num__covered  ■*  1? 
for  (1-0;  i  <  nvar;  ill) 

num_covered  *-  (I->B f i ). upper  -  I->B ( i | . lower  4  1); 
if  ( (num_xeroed  >  cur_num_zetoed)  (f  ( 

(num  zeroed  — —  cur_num__*etoed)  (fi 
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n 


i 

i 

i 

i 

return (SI_best) 


(num  covered  »  cur_num  covered) 

( 

curntimreroed  —  nun  retoed; 
cur_num_covered  -  num  covered; 
copy^implicant (SI_best,  I) ; 
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/*  comnon ,  c 


-  Heuristic  support  functions  common  to  DtM  and  F*A. 

*/ 

•include  "defs.h" 

int  *eval  (E, x) 

Expression  *E ; 
int  x; 

/*  - 

: function: 

-  Evaluate  the  expression  at  (xl ,  x2,  x3,  x4 ,  x5 .  ,  .  . xn> 

-  Builds  a  vector  of  coordinates  for  _eval(( 
scalls: 

_eval  (1 
: returns : 

-  A  vector  containing  the  Integer  result  0  <—  eval  ()  <~  radix 
and  a  flag  set  if  the  value  reached  the  highest  logic  value 

( 

int  *X; 

X  -  (int  *)  <Sx,‘ 
return (_eval (E, X) I ; 

I 

Int  *_eval (E, X) 

Expression  *E; 
int  *X; 

/* . . - . - . . . . . . . 

: function : 

-  Evaluate  the  expression  at  X,  where  X  is  a  vector  of  coordinates 
iglobals : 

nvar 

radix 

! side_ef  fects : 

STAT 

{called  by: 

mlmO  -  pa.c 

valid^lmpl leant  0  -pa.c 
plck__lmpllcant  0  -  pa.c 
mim<)  -  dm.c 
valid  ImpllcantO  -  dm.c 
pick  ImpllcantO  -  dm.c 
_cfO 

comrute_rbc  0 
eval  () 

gen_bounds  0 
prlnt_map  () 

: returns : 


-  A  vector  with  the  value  of  the  expression  at  the  specified 
coordinate  as  its  first  element,  and  a  flag  set  if  this  value 
has  attained  the  highest  logic  value  (IlLV) 


int  nterm  -  E->nterm; 
register  i,  j,  k; 
int  out_of_bounds; 
static  int  V(2|; 
register  rml  —  radix-1,* 

STAT->calls_eval++; 

V ( EVALJ  -  0; 

V(HLV)  -  0; 

/*  for  each  term  ...  */ 
for  <i«0;  1  <  nterm;  ltt)  | 

/*  for  each  variable  ...  */ 

for  ( j-0, out_of_bounds-0;  j  <  nvar;  j++)  ( 

if  ( 

(X(j|  <  E->T(i|.B(i|. lower)  II 
(*ljl  >  E->TU1  .Blj)  .upper) 

)  I 

out_of_bounds  ■*  1; 
break; 

I 

I 

if  (out_of_bounds) 
continue; 

/*  if  this  is  a  don't  care,  return  the  radix  */ 
if  (E->T(i) .coeff  —  radix)  ( 

V(EVAL)  -  radix; 
return (V) ; 

} 

V[EVAL|  +-  E->T(i) .coeff; 
if  (V)EVALl  >-  rml)  ( 

/*  set  a  flag  which  means  E_orlg  was  saturated  at  this  X  */ 
V  |  HLV  J  -  1; 

I 

if  (V(EVAL)  >  rml)  | 

V ( EVAL )  -  rml; 

I 

else  If  (VIHLV)  SS  <V{EVAL)  <-  0()  I 
V | EVAL)  -  radix; 
return (V) ; 

I 

) 

return (V) ; 
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int  *next_coord (T, first) 

Term  *T; 
lnt  first; 

/* - - - - - 

: f unct ion : 

-  Compute  the  next  possible  coordinate  for  term  *T 

-  If  first  — 1,  initialize  the  coord  vector 
:called_by : 

mim()  -  pa.c 
va 1 idlmpl leant ( )  -  pa.c 
mim  ( )  -  dm. c 

valid_impllcant ()  -  dm.c 
compute_rbc  0 
: returns ; 

-  An  integer  vector  containing  the  coordinates. 


( 

static  int  coord (MAX_VAA4 1 ) ; 
static  i; 

/*  if  the  first  time  through,  load  the  vector  */ 
if  (first)  ( 

for  (i-0;  i  <  nvar;  144)  | 

coord(i)  —  T->B ( i J . lower ; 

J 

) 

else  ( 

i  -  0; 
coord ( 1)44; 
for  <;;)  ( 

if  (coord[i|  >  T~>B|1 J .upper)  ( 
coord (ij  -  T~>B|i) .lower, 

144: 

if  (i  >-  nvar) 

return (NOLL) ; 
coord (i) 44; 

I 

else  ( 

break; 

) 

I 

1 

return (coord) ; 


I 


Pound  *gen_bounds (X) 
int  *  X ; 

/* - 

: f unct ion : 

-  Generate  the  permissible  bounds  around  location  X  in  the 
working  expression 

: global s : 
radix 
nvar 
E_work 
E_or ig 

: side_ef  f ect s : 

STAT 

: cal led_by : 

pick_implicant  ()  -  pa  .c 
pick_impli«jant  ( )  -  dm.c 
: calls : 

_eval 0 
vcopy  0 
: returns : 

-  A  bounds  array 


Static  Bound  B  (tlAX_VAR+ 1  J  ; 
int  nterm  —  E_work . nterm; 
registe  i ,  j , k ; 
int  value,  Vw  1 2) ,  Vo 1 2) ; 
int  XptHAX_VAR+lJ ; 

value  —  X(nvar); 

STAT->cal Isgenbounds 1 4 ? 

/ *  for  each  variable  (direction) ...  V 
for  (i-0:  i  <  nvar;  i  +  4  )  ( 

/*  dup  the  coordinate  */ 

for  (3-0?  j  <  nvar;  j  +  D  Xp(j)  -  Xtj); 

Bli) . lower  -  X[iJ ; 

/*  while  not  on  a  left  hand  edge,  move  left  V 
while  (Xp(iJ  >  0)  ( 

Xp(i|— ; 

vcopy (Vw,_eval (5E  work, Xp)  )  ? 
vcopy (Vo, _eval ( ftF_or ig, Xp) ) ? 

/*  if  can't  expand  to  left  ....  */ 

if  ((((value  >  Vw  ( EVAL  J )  t£  (VolEVAL)  <  (radix-1)  )))  ( 

B(i|. lower  -  Xp ( i ] ; 

) 

else 

break ; 

) 
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/*  dup  the  coordinate  */ 

for  (j-0;  j  <  nvar ;  j++)  Xp(j]  -  XI j); 

B  [i J  , upper  -  X l i J ; 

/*  while  net  on  a  right  hand  edge,  move  right  */ 
while  (Xp(i)  <  (radix-1))  ( 

Xp [ i ) +  + ; 

vcopy (Vw, _eval (&E_work, Xp) ) ; 
vcopy (Vo, eval ( &E_orig,  Xp)  )  ; 

/*  if  can't  expand  to  right  ...  */ 

if  ((((value  >  Vw(EVAL))  &&  (Vo(EVAL)  <  (radix-1))))  ( 

B l i ] .upper  -  Xp(i) ; 

) 

else 

break ; 


I 

return  (B) ; 
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/*  Working  structures  for  picking  the  next  implicant  within  bounds  * / 


static  Bound  I B (MAX_VAR4l ) ;  /*  Current  bounds  •/ 

static  Term  I;  /*  Implicant  */ 

static  int 
I_var, 

Ifirst, 

Ival ; 

int  X_origlMAX_VAR+l J ;  /*  Where  we  start  */ 


init_impl leant (X) 
int  *X; 

/* - - - - - 

: f unct ion : 

-  Initialize  the  static  term  structure  above  from  which  successive 
implicants  will  be  returned 

-  X  is  the  starting  minterm 
: slde_ef  f ect s : 

-  The  structures  above 
! called^by : 

pick_imp] leant  0  -  pa.c 
pick_implicant ()  -  dm.c 


int  nterm  -  E_work . nt erm; 
register  i; 


/*  initialize  the  Implicant  */ 

I.B  -  IB; 

I.coeff  -  X  Invar); 

for  (i«^0;  i  <  nvar;  i.  4  4 )  ( 

Z . B ( i 9 .upper  -  Xfi J; 

I.B til .lower  -  X(iJ ; 

) 

Ivar  -  0; 

Ifirst  -  1? 
t_val  —  x I nvar ) ; 

for  (i-0;  i  <  nvar;  1+4)  X_origliJ  -  X(i); 


Term  *next_lmplicant (B) 

Bound  *B; 

/*  - 

: f unct ion : 

-  On  each  call,  return  the  next  implicant  within  bounds  B 
: side_ef f ects : 

STAT 

: called_by : 

pick  impl i cant  ()  -  pa.c 
pick_implicant ()  -  dm . c 
: returns : 

-  An  implicant  as  a  term  structure 


Int  nterm  -  E_work . nterm; 
int  Xp(MAX_VAR  +  l J  ; 

9TAT->cal 1 s_next_impl leant ++; 
if  flMrst)  I 
I_flrst  —  0; 
return  (SI) ; 

I 

while  (I_var  <  nvar)  ( 
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/*  expand  left  */ 

X  .  T ;  I_  var)  .lower--; 

/*  If  we  can't  go  further  left,  then  ...  »/ 
if  (I . B ( l_var 1 . lower  <  B | I_var) . lower)  ( 

/*  move  back  and  go  right  */ 

I . B ( l_var | . lower  -  X_orig | l_var ) ; 

I . B | I_var  J . upper +  +  ; 

/*  if  we  can't  go  further  right,  then  ...  */ 
if  (I . B ( I_var ] . upper  >  B ( I_var ). upper )  ( 

/*  reset  and  go  to  the  next  higher  dimension  •/ 
I . B ( I_var) . upper  -  X_or ig I I_var ] ; 

I_va  r 1 4  ; 
continue; 

I 

I 

l_va  r  —  0 ; 
return (SI) ; 

) 

return  (NULL)  ; 
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int  copy  impl leant (dest, ace) 

Term  *dest,*9rc; 

/*  - 

: f unct ion : 

-  Copy  the  implicant  pointed  to  by  etc  to  dest 
: called_by : 

plck_impl leant ()  -  pa . c 
pick  Implicant ()  -  dm.c 


register  i; 

deat->coeff  —  arc->coeff; 
for  (i-0;  1  <  nvar;  i44)  ( 

dest->B  {  i  |  .  lower  —  src->B  f  i  ).  lower  ; 
deat->B ( i | . upper  -  src->B ( i ). upper ; 

I 

I 


subtract_lmpllcant (1) 

Term  *  I ; 

/*  - - - 

: function : 

-  Add  implicant  I  to  the  working  expression  as  a  negative  term 
(negated  coefficient) 

-  Add  implicant  1  to  tthe  final  expression 
: global s ! 

HEUR 

nvar 

: side_ef fects : 

Ework 
E  final (| 


register  l.term; 

term  -  E_work . nterm; 

E_work . nt erm4  4 ; 

Ework.T  -  alloc_term (Ework . T, - (I->coeff) , E_work . nterm) ; 
for  (1-0;  i  <  nvar;  144)  | 

E_work .T(term| . B { 1 J . lower  »  I->B ( 1 J . lower; 

E_work.T(termJ  . B | i. ]  .upper  ■»  I->B ( i ). upper ; 

I 

term  -  E_f  inal  (HEUR)  .nterm; 

E_f  Inal  (IIEtlRJ  .nterm+4  ; 

E_f  Inal  (HEUR)  .T  -  alloc_term(E_f  inal  (IIEUR)  .T,  I->coef  f,  E_f  inal  (HEUR)  -  nterm)  ; 
for  (1—0;  i  <  nvar;  144)  | 

E_f Inal  I HEUR ) . T l term) . B ( 1 )  . lower  —  I->B(l).lower; 

E_f inal (HEUR) .T (term) . B [i ). upper  -  I->B ( i ). upper ; 

I 

I 

/*  vcopyO 

-  copies  the  value  vector  from  a  to  d 

»/ 

vcopy (d, s) 

Int  *d, *s; 

( 

d l 0 |  -  9(0); 

dll!  -  »)l); 

I 
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•Include  "defs.h 


/ *  memory  allocation  functions - - - - - - - - - */ 

Term  * a  1 locate rm (p, coe f f , n) 

Term  *p; 
lnt  coeff,n; 

/* - - - - - 

! function : 

-  Allocate  apace  for  a  term  array,  initializing  the  last  element 

-  if  p  la  HULL,  allocate  new  apace 

-  If  p  la  not,  realloc 
j  returns : 

-  A  pointer  to  the  Term 


char  *mallocO,  *realloc(); 

Bound  *alloc_bound ()  ; 

if  (p  —  NULL)  ( 

if  Hp*(Tetm  *  )  mal  loc  ( si  zeof  (Term)  4n)  )  ***•  NULL) 
fatal ("alloc_term {) :  out  of  memory\n") ; 
p->coeff  ~  coeff? 
p->B  alloc_bound  { )  ? 

I 

else  | 

if  ((p-|Term  *)  realloc (p,  si zeof (Term) *n) )  —  NULL) 
fatal ( "alloc_Tetm () j  out  of  memory\n") ; 
pln-1). coeff  -  coeff; 
p(n-l).B  -  alloc_bound  ()  ; 

I 

return (p) ; 

1 

Bound  *alloc_bound 0 

/* - - - 

?f unction: 

-  Allocate  space  for  F.^orig.nvar  bounds  entries  and  initialize 
each  bound  to  -1,E  or ig . radix-1 . 

-  If  p  is  NULL,  allocate  ne*  space 
: global s : 

F_or ig 
: returns : 

-  A  pointer  to  the  Bound  array 

- - -  v 

1 

Bound  *p; 

char  *malloc(); 
register  i; 
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if  (  (p—  (Bound  M na 1 loc ( ai leof ( Bound) * (E  orig.nvsr) 
fatal ("alloc_bound  O !  out  of  nemory\n") ; 

for  (i-'-O;  i  <  E_orig .  nvar ?  i++)  ( 

p | 1 J , low?r  •  -1 ; 
p ! 1 1 . upper  -  E_or ig . radix- 1 ? 


return  (p) ; 


f  >  —  null; 
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init_expr  () 

/*  - 

: function : 

-  Initialize  E_work,  E_orig  end  E_final 
: side_ef  f ects : 

Ework 
E_orig 
E_f inal 

- - -  ./ 

I 

Evork.T  -  HULL; 

E_orig.T  -  NULL; 

E_orlg.nvar  —  0; 

Eorig.nterm  —  0; 

Eorig. radix  -  0; 

E_f Inal (0) .T  -  NULL; 

E_f inal ( i ) . T  -  NULL; 

I 


deail  oc_expr (e) 

Expression  *e; 

/«  - 

: function: 

-  Deallocate  the  expression  pointed  to  by  e 


( 

Term  *p; 
register  i; 

if  <e->T  1-  HULL)  | 

for  (p  -  e->T, i-0;  i  <  e->nterm;  i++) 
if  <p C 1 1 . B  I-  HULL) 
free  (pH)  .B)  ; 

tree (p) ; 

) 

I 
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dup_expr (E_dest , E_src) 

Expression  *E_dest, *E_src? 

/* - 

: function : 

-  Duplicate  the  expression  pointed  to  by  Esre  by  allocating 
and  copying  into  the  expression  pointed  to  by  E_dest . 
scalls : 

alloc_bound  ( ) 

-  */ 

( 

Term  *T; 

Bound  *B; 
register  i, j; 
char  *malloc(); 

E_dest->radix  -  E_src->radix; 

E_dest->nvar  -  E  src->nvar; 

E_dest->nterm  -  E_sr c->nterm; 

If  ((T-(Term  * ) malloc (sizeof (Term) * (E_dest->nterm) ) )  —  HULL) 
fatal < "dup_expr 0 :  out  of  memory\n"); 

for  (i-'O;  i  <  E_src->nterm;  !♦  +  )  ( 

T|i).coeff  -  E_src->Tf i) .coef f ; 

T ( i ) . B  -  alloc_bound ( ) : 

for  (j-'O;  j  <  E_src->nvar;  jit)  f 

Tf i ) . B | j) . lower  -  Esr c->T f JJ . B | j) . lower ; 

T(i| .  B( j| .upper  -  E_src->T l i | . B [ j ) . upper ; 

» 

) 

E_dest->T  -  T; 

I 
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/*  mvltest.c 


-  A  teat  source  file  generator 

*/ 


llnclude  <stdio.h> 


chnr 

•version 
•usage  “ 

•help tJ  ■ 


HULL 

If 


-  "NntIVLTEST  Vl.O  Copyright  1988  Naval  Postgraduate 
"usage:  mvltest  (-SN  -rN  -IN  -tH  -oN] \n", 

■  I 


-sH  -  Cycle  the  random  seed  H  times  (default 
—  r H  -  Set  radix  to  H  (default  «•  4) ", 

-IN  -  Set  number  of  Inputs  to  H  (default  *■  2) 
-tN  -  Set  number  of  terms  to  H  (default  “  1)" 
-oH  -  Set  number  of  output  to  H  (default  —  1) 


0)  ' 


School \n\n" 


main (argc, argv) 
char  *argv|f; 
( 


Int  nterms  *  1; 

Int  rers  -  1, 
radix  —  4, 
nvar  —  2, 
fcoeff  -  1; 
int  seed  0; 
int  i, j, k; 

int  *b,  *get_bound  ()  ; 
float  random(); 
char  *p; 
register  ap; 


for  (ap~l;  ap  <  argc;  ap++)  ( 

p  -  argv |ap| ; 
if  <*P14  -  '-•)  ( 

while  <*p)  ( 

switch  (*p++)  ( 

case  ' r’  : 

sscanf (p,  "td", (radix) ; 

*P  -  ' \0' ; 

break; 
case  ' 1’  : 

sscanf (p,  "%d", (nvar) ; 

*P  -  '\0'  ; 

break; 
case  ' t' : 

sscanf (p,  "%d", (nterms) ; 
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I 

for 

for 


*P  -  ' \0' ; 
break; 
case  ' o' : 

sscanf (p, "Id",  Steps)  ; 

‘p  -  ' \0- ; 
break ; 
case  ’ s'  : 

sscanf (p,  "Id", Sseed) ; 

*P  -  '\0'  ; 

break; 
case  ' : 

fprlntf (stderr,"lsls", version,  usage)  ; 
for  (ap-O;  helplap);  ap4  4) 

fprlntf (stderr,"ls\n",help(ap)) ; 
exit (0) ; 
default : 

fprlntf (stderr, "mvltest :  unknown  switch  -  %c\n”,p[-l)) 
fprlntf (stderr, "Is", usage) ; 
exit  (1)  ; 

t 

) 

I 

(1-0;  1  <  seed;  144)  random(); 

(1  -  0;  1  <  reps;  144)  | 

prlntf ("%d:ld: ",  radix, nvar) ; 
for  (j-0;  )  <  nterms;  }4  4)  | 

prlntf ("\n  4ld",  (lnt!  (random!) * (float)  (radix-1) ) 41) ; 

for  (k-1 ;  k  <-  nvar;  k44)  ( 

b  -  get  boond(radix); 
prlntf ("‘Xld (Id, Id) »,k,b(0) ,b(l) ) ; 

I 


I 


prlntf (";\n") ; 


lnt  *get_bound (radix) 

( 

static  lnt  b(3); 
float  random(); 

for  (;;)  ( 

b(0)  —  (lnt)  (random!) ‘((float)  (radix) -0 . 001 ))  ; 
b(l)  -  (lnt)  (random()  *(  (float)  (radix)  -0.001)  )  ; 
if  (b ( 0)  <-  b(l)» 
break; 

) 

return (b) ; 

I 

float  random!) 

( 

static  long  a-100001; 

a  -  (a‘125)  %  2796203; 
return ( (float) a/2796203) ; 

I 
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I  Makefile  for  mvl 

I 

erLAGS— o 

OBJS--main.o  dm . o  pa.o  common. o  parse. o  alloc. o 
all :  mvl  mvltest 
mvl:  $ (OBJS) 

cc  $  (CfLAGS)  $ (OB JS)  -o  mvl  -11 

mvltest :  mvltest.o 

cd  $  (CFLAGS)  mvltest.o  -o  mvltest 

ma 1 n . o:  ms in.c 
alloc. o!  alloc. c 
parse. o!  parse. e 
dm . o :  dm . c 
pa.oi  pa.c 
common . o !  common. c 
parse. c!  parse. 1 
lex  parse. 1 

mv  -f  lex.yy.c  parse. c 
mvltest.o:  mvltest. c 
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